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REPLACEMENT WHEN CONSTANT 
FAILURE RATE PRECEDES WEAROUT 


Larry Hunter and Frank Proschan* 


General Telephone and Electronics Laboratories, Inc., Menlo Park, Calif. 
and 
Sylvania Electric Products, Inc., Mt. View, Calif. 





Assume a component has a constant failure rate during the first T unit 
of time following installation. After that, wearout sets in so that fail- 
ure rate begins to increase sharply. Assume further that replacement 
is made at the time of component failure or at T units of time following 
installation, whichever occurs first. Under these assumptions, we 
obtain explicit formulas for the distribution and expected value of the 
number of planned replacements, the number of failures, and the total 
number of removals due to either planned replacement or failure 
replacement, corresponding to any specified length of time. Examples 
are worked to show the application of these formulas. The results 
obtained in this article are of use in determining the number of spares 
to stock or the budget required to maintain the equipment. 











INTRODUCTION AND SUMMARY 

In many reliability problems, a reasonable assumption is that component failure rate is 
constant or equivalently, that component life follows an exponential distribution, until wearout 
begins. A procedure commonly followed is to schedule a replacement at that time so as to 
reduce the probability of component failure during actual operation. In this situation it is of 
value to know for any specified length of time, the distribution and the expected value of the 
number of planned replacements, the number of failures, and the total number of removals due 
to either planned replacement or failure replacement. This kind of information is needed, for 
example, in determining the number of spares to stock or the budget required to maintain the 
equipment. 

In this article, we derive these distributions and expected values, and give examples to 
illustrate the results. 


PRECISE STATEMENT OF THE MODEL . 
We assume component life is governed by the density g(t) = e' for 0 St S7 (or 


t 
equivalently the failure rate (t)/1- ®(t) is constant, namely 1, where &(t) = f ¢~(u) du). 
0 


After time T the component begins to wear out so that the failure rate is greater than 1; it will 
not be necessary to specify the exact form of the distribution for t > 7 since the component 
will be immediately replaced upon failure or after 7 units of time, whichever occurs first. 
There is no loss of generality in assuming the parameter of the exponential to be 1; essentially, 
this amounts to measuring time in units of exponential mean life. 


*Now at Boeing Scientific Research Laboratories, Seattle, Washington 
4 


















L. HUNTER AND F. PROSCHAN 






Starting with a new component at time 0, let 
Ep (t) = expected number of removals during [0,t] due to either failures or planned 
replacements (""removals" will refer hereafter to both planned and failure 
replacements), 


R(n,t) = probability of at least n removals during [0,t], 

Ep (t) = expected number of planned replacements during [0,t], 
P(n,t) = probability of at least n planned replacements during [0,t], 
Ey (t) = expected number of failures during [0,t], 

F(n,t) = probability of at least n failures during [0,t]. 


Next, we derive closed expressions for these quantities. 


NUMBER OF REMOVALS 


The process of making removals is clearly a renewal process. By standard renewal 
theory arguments ([1], p. 252): 














* 

(3.1) ER (s) = =, 
* 
1-r (s) 

where 

fe @) 

Ep (s) = f et dEL(), 

and 





. 
r* (s) = [ e7st y(t) dt + e“*11_6(7} . 





- 
(Note that r (s) is the Laplace-Stieltjes transform corresponding to the probability distribu- 
tion for removal by time t.) Since ¢(t) = e~* and 1- (7) = a we obtain 





oe -T (s+1) 
r(s) = {1+ 50 }. 


Thus from (3.1), 





Ep (s) = {2 +e tt {1 -e a sa ° 


Inverting the Laplace transform using [2]: 





(3.2) 






Ep (t) =t 
O=nSt/t 


FAILURE REPLACEMENT 


0 when t < (n+1)7T 
1 when t 2 (n+1)T 


Equivalently, for (k- 1)7 St <k7, we have 


1 


(1 - a 


(3.3) Ep, (t) = 


{t - @+r-e™* - eT . (t+1-k7)e7*7 + (t+1-(k- 1) 7yeket)r| 


for k= 1, 2,.... See Figure 1 for a plot of Ep (t) for T= 1 and T= 2, 





























Figure 1 - Expected number of removals in time t 


Thus (3.2) or (3.3) gives a closed form for the expected total number of removals 
during [0,t]. 
Note that for small s, 


* ami=k << 
Ep(s) ~{1 - e7’} s*. 
Hence by the Tauberian theorem of ([3], p. 192) 


(3.4) Ep(t) ~ {1 - e Vie t 


for large t. (See also Ref. [1], p. 246.) This is readily confirmed from (3.3). 
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To obtain a closed form for R(n,t), the probability of at least n removals during [0,t], 
we use the fact that the nt th power of the Laplace-Stieltjes transform corresponds to an n-fold 
convolution of the random variable: 





© n -jt(s+1) 
-st Pe n j healed 
! e dR(n,t) = r (s) “2, (;) s ene . | 
Thus by Ref. [2], if (k-1)7 =t << kT 
t n-1 min (k-1,n) j-1 - st 
4 n\ a (t-jr) 
-6 -t 
R(n, t) = | e (n+1)! dé + a (oa { (n-1)! . 





—_ 
" 
_ 





Using the Euler formula for the repeated derivative of a product, we get 





t n-1 min (k-1,n) j-1 
_{ 0 8 M\ +t - ) (t-nn) 
(3.5) R(n,t)= | e cae+ 2 (;)e* D> eam (n-j+m)! 
0 j=l m=0 . 


for (k-1)7 =t <kt, k=1, 2,.... Note that the integral in (3.5), representing the incomplete 
gamma function, may be alternately expressed as 


>. P 
- e~? 4) 








and obtained numerically from the tables in Ref. [4]. See Figure 2 for a plot of R(n,t) for 
T=1 and T= 2, 


NUMBER OF FAILURES 
F(n,t), the probability of at least n failures during [0,t] and Ep (t), the expected 
number of failures during [0,t], are immediately obtained by noting that the probability of 
failure during any interval of length dt is equal to dt + o(dt), and is independent of the number 
of preceding failures or of the age of the component currently operating. Hence 





n n+1 
“<< /: % 

4.1 rah a€" — «+ ee? ——— 

ate (n, t) n! (n+1)! 


which is the right-hand tail area of the Poisson distribution ([5], pp. 401-402); F(n,t) is 
tabulated in Table II of Ref. [6]. 
We obtain immediately 


(4.2) 





Ep (t) =t, 


the mean of the Poisson distribution. 
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Figure 2 - Probability of at least n removals in time t 
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NUMBER OF PLANNED REPLACEMENTS 
To obtain Ep (t), the expected number of planned replacements during [0,t], we simply 
note that 


Ep (t) = E,(t) + Ep (t) . 


Hence, we immediately obtain 


1 —— . 
Ep(t) = ~Aazite (t+ 7-1) e°7 - e727 - (t+1-k7) eo” 


i-@” 


(5.1) + (t+1-(k- 1) 7) e+) 7} -t, 


using Eqs. (3.3) and (4.2). 
To obtain P(n,t), the probability of at least n planned replacements during [0,t], let 
us define 


fo @) 
Ep(s) = j eS GE,(), 


A(t) = probability of planned replacement by time t , 


’ 


fe @) 
A*(s) = I e~St a A(t). 


Again using the standard renewal theory result ([1], p. 252), 


A’ (s) 
* Ss 
Ep(s) = ——,—, 
1-A (s) 
we have 
* 
Ep (s) 
hese 
1+ Ep (s) 
But we know 
c 





Ep(s) = Ep (s) - l t eS at 








- {2 ‘ ideas {i 7 eT sean} - 3 s *s+1 


s se 7(Stl) “a . 
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nply 
Thus, we have immediately 
= n 
- * (s + 1) 
P’ (n,s) = | eSt gP(n,t) = A (s)" = = 
, {se 7(s+2) + i} 
To invert P” (n,s), we write 
a se 
P*(n,s) = e7NT ». (5) gin? _-nTs 
let j=0 


n 
- ne (ntl) 7 . ginn-1 .-(n+1)7s 
AG 


n 
n(n+1)  -(n+2)7 M\ j-n-2 _-(n+2)7s 
+ ~s e 2, (|) s) e + see 


Then we obtain from Ref. [2] 


n 
= emt (t - nz)" 
P(n,t) = e” 2; fea 





ne7 (n+l) T a % (t - (n+ 1) or j 





= (n + 1 - j)! 
n 3-4 
n(n+1) _-(n+2)7 N\ (t - (n+ 2) 7)"*OI 
an ——— = Zi) inch ii tees 


j=0 
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Equivalently, for (k- 1)7 =t <k7T, we have the closed form 


lV 


k 


(0 for n 


n 


_ a an-j ; a z n+1-j 
eat » bh — ~ new (ntl) Tt a (3) (t ae 





j=0 j=0 





n . 2-j 
n(n +1) _-(n+2) DY (t - (n+ 2) 7/2) 
(5.3) P(n,t) =< * ~ oy : , 2, (5) (n + 2 - j)! 


+ (eykeden n(n +1)... (k- 2) ,-(k-1)7 
(k- 1-n)! 





My (2) (t= (k - 1) 7) 
oF be gee ae forn<k, 


[ z 


where k= 1, 2,... . See Figure 3 for a plot of P(n,t) for 7=1 and 7 = 2, 

To adapt any of the results (Eqs. (3.2), (3.3), (3.4), (3.5), (4.1), (4.2), (5.1), (5.2), or 
(5.3)) for the general case (t) = re" , 0 =t=T, for component life, simply replace t by tA 
and T by TA. 








EXAMPLE 

A tunable magnetron in an electronic countermeasures transmitter has a constant fail- 
ure rate of A = 1/10 per hour for the first 10 hours of operation. After 10 hours of operation, 
the rate of failure increases rapidly due to wearout of the mechanical tuning components. To 
decrease the probability of failure during actual battle, replacement of the magnetron is 
scheduled after it has been in operation for t' = 10 hours. To give management a complete 
picture of spares requirements and possible operating difficulties, determine the expected 
number of actual failures, planned replacements, and removals due to either failure or planned 
replacement during the first t' = 35 hours of operation. In addition, determine the probability 
of at least n failures, of at least n planned replacements, and of at least n removals due to 
either failure or planned replacement, during the first t' = 35 hours, where n= 1, 2, 3,.... 

Note first that 7'A = 1 and t'A = 3.5. Hence, to obtain the expected number of actual 
failures during the first 35 hours, we use Eq. (4.2) with argument t = 3.5 to find E , (3.5) = 3.5. 
To obtain the expected number of removals (due to failure or planned replacement during the 
first 35 hours), we read off an ordinate of 5.2 corresponding to an abscissa of 3.5 on the curve 
labelled 7 = 1 in Figure 1 (based on Eq. (3.2) or Eq. (3.3)). To obtain the expected number of 
planned replacements during the first 35 hours, we recall that Ep (3.5) = Ep (3.5) - E 7 (3.5) 
= §.2 - 3.5 = 1.7. 

Next, to obtain the probability of at least n removals in the first 35 hours (due to either 
failure or planned replacement), we use the curve labeled 7 = 1 in Figure 2, based on Eq. (3.5). 
For example, for n= 1, 2, or 3 we read off a probability of 1; for n = 4, we read off a proba- 
bility of .9. To obtain the probability of at least n planned replacements during the first 
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Figure 3 - Probability of at least n planned replacements in time t 
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35 hours, we use the curve labeled 7 = 1 in Figure 3, based on Eq. (5.2) or Eq. (5.3). For 
example, P(1, 3.5) = .94, P(2,3.5) = .62, and P(3,3.5) = .15. Finally, to obtain the probability 
of at least n failures during the first 35 hours, we use Eq. (4.1) which is tabulated in Table I 
of Ref. [6]. Thus F(1,3.5) = .97, F(2,3.5) = .86, F(3,3.5) = .68, F(4, 3.5) = .46, etc. 
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THE ORGANIZATION AND OPERATION OF A TAXI FLEET 


Richard F. Meyer and Harry B. Wolfe 


Arthur D. Little, Inc. 
Cambridge, Massachusetts 





In dispatching a taxi fleet to meet random customer demand, two types 
of strategies must be defined: (a) the dispatching strategy, which sets 
the rules for choosing a taxi to serve a given call, and (b) the idle time 
strategy, which determines how the idle taxis are to be distributed 
most advantageously to meet future demand. The effect of dispatching 
and idle-time strategies on taxi utilization and on passenger service is 
demonstrated under several conditions of operation. First, the simplest 
taxi model is examined: random demand between two points serviced by 
one taxi. Next, we consider the effect of different strategies on the 
operation of a taxi fleet when demand is uniformly distributed in a two- 
dimensional region. In the last section, the main properties of these 
simplified models are exhibited in an intermediate case by a Monte 
Carlo simulation of taxi service between a number of discrete points. 
The purpose of the paper is to provide simple and practical methods to 
determine the payoff in customer waiting time and/or taxi savings 
resulting from more sophisticated dispatching and idle-time strategies. 
We conclude that good taxi redistribution policies markedly decrease 
customer waiting times, except when the level of taxi utilization is too 
high. Our results apply most directly to private transportation fleets, 
rather than to automobile taxi companies where cruising may be 
important. 











INTRODUCTION 

Based on an actual case history, this article discusses the scheduling problems arising 
in the design and operation of a transportation system. The actual system provides transporta- 
tion for varied demands, personnel and material, scheduled as well as random, between points 
of a two-dimensional region. The dominant type of demand, however, is for unscheduled trans- 
portation; therefore, we discuss only the problems arising in the design and operation of what 
might be called a taxi fleet. 

The case being considered was not subject to the restrictions of streets; the "taxis" 
could always travel fairly directly between any two points. Also there was no possibility of 
picking up passengers by cruising. In these respects, the operation differs from that of the 
usual automobile taxi fleet and corresponds, for example, to a fleet of helicopters which are 
on call to transport personnel and material between points of a large two-dimensional region. 
The objective of this helicopter taxi system is to provide a predetermined level of passenger 
service, with the least number of helicopters. In order to apply to automotive taxi fleets, our 
analysis could be extended to include picking up passengers by cruising. 
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The policies that may be considered in the dispatching of this helicopter taxi fleet are: 

(a) When a call arises for transportation, which taxi should be assigned to service the 
call? This we call the dispatching strategy. 

(b) When a taxi is free, where should it be relocated most advantageously with refer- 
ence to future demand? This is the idle-time strategy. 

We wish to determine the following: How do the dispatching and idle-time policies affect 

passenger service time and taxi utilization? How are these affected by the distribution of 

demand? How does passenger service depend on the number of taxis? 

In the next section, we investigate the effect of different scheduling policies for the 
simplest possible model: random demand between two fixed points to be served by one taxi. 
This model is simple enough to be treated analytically and provides considerable insight into 
the situation. When the number of points of demand and the number of taxis increases, it 
becomes difficult to treat the problem analytically. However, if the number of points of 
demand is very large and reasonably uniformly distributed, the distribution of demand may be 
approximated by a uniform distribution density, and analytical treatment again becomes feasi- 
ble. This article also discusses the effect of different strategies on the operation of a taxi 
fleet which is serving uniformly distributed random demand in a large, two-dimensional region. 
Later, we obtain similar results for an intermediate case by a simple Monte Carlo simulation 
of taxi service between several demand points. 


MODEL OF TWO SOURCES SERVED BY ONE TAXI 

The effects of different idle-time strategies can be seen from studying a simple model 
in which two sources of taxi demand are served by one taxi which moves back and forth between 
them. Ferry-boat operation, shuttle service between two plants, helicopter service between 
two airports, etc., can all be approximated by this model. 

The two sources are shown in Figure 1, located respectively at Points 1 and 2 with 
average demands of aT and Vo calls per hour, and a travel time 7 between them. A taxi which 
is occupied by a passenger, or is traveling (empty) in order to pick up a passenger, will never 
be allowed to turn around to pick up another passenger until it has reached its destination. Two 
simple idle-time policies will be compared: 

Policy A. A free taxi always remains at the place where it becomes free. 

Policy B. A free taxi always returns to the busier of the two points. 


"2 "2 
vy, V5 Figure 1 - Shuttle service, showing ratio 
between Points 1 and 2 


oa. bs — 
oe] oOo 





On comparison, it will be found that Policy B is generally the better one from the point of view 
of passenger waiting time, except when the taxi is operating at a very high level of utilization, 
in which case it is found that the additional "deadheading" (traveling empty) required under 
Policy B is uneconomic in terms of total passenger waiting time. 

Policy A. A free taxi remains where it becomes free. Assuming that the calls arising 
at Points 1 and 2 occur randomly in time and independently of one another, the probability of 
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having to make two trips from Point 1 to Point 2 in succession, without an intervening trip 
from Point 2 to Point 1, is: 


co io) vi? 


-Vt' . 
(1) { v, dt { vie dt’ = where V=V,+V5. 
, 8 ss y2 -”*S 





The corresponding probability of two trips from Point 2 to Point 1 in succession is 
similarly v9"/ y2 so that the probability of having to make two similar trips of either type in 
a row is 


2 
Vv + V 
1 2 
(2) ——__—_—.. 


y2 


If two similar trips are made in a row, the taxi deadheads back to pick up the second passenger, 
and the second trip therefore requires a running time 27 instead of tT. The remainder of the 
time, i.e., a fraction 





of all trips, the taxi gets alternating calls, and therefore has no need for deadheading. If F is 
the total running time per unit time, then: 








Yo” RY 2vav Viv 
1 
(3) P 4 tency Bet -* St Sr 2) or, 
y2 v2 y? 
For simplicity, using the notation 
p= v,/v q=1- p= v/v X= VT, 


we may rewrite Eq. (3) more simply as 
(4) F = 2(1 - pq)x. 


x is the fraction of time that passengers are being transported and, hence, is called the 
"utilization." 
The average waiting time for a passenger from Point 1 consists of two parts, as follows: 
(a) If the taxi is stationary (probability 1-F), the passenger must wait a time 7 if the 
taxi is at Point 2 (probability v,/ v), and will not wait if the taxi is at Point 1; hence, 
this contributes (1-F) (v,/ Vv) T to his expected waiting time. 
(b) If the taxi is running (probability F), the expected waiting time is 7, since the taxi 
makes the same number of trips in either direction; hence, the contribution is FT. 
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A passenger from Point 1, therefore, on the average waits a time 


< 
wn = T + r*| ‘ 


A similar expression holds for a passenger from Point 2. The total passenger waiting time 
under policy A, T A? is now found by adding the total waiting time for passengers from Point 1 
to the total waiting time for passengers from Point 2: 


= Vy an "1 T + P| + Vo on “ + F 
v 


> 
i} 


Vv 


(5) 


i} 


(1 - 2pq)x + 4pq(1 - pq) x”. 


In this approximation, the total passenger waiting time is a quadratic function of the 
level of utilization x. 

Policy B. A free taxi returns to the busier point. Arbitrarily, consider Point 1 as the 
busier point, so that p is larger than q (v, > Vo). The detailed rule for the idle-time strategy 
then becomes: as soon as a taxi becomes free, it must always return to Point 1. We call this 
procedure "idle deadheading."' There are two simple idle-time policies for handling calls 
which arise at Point 2 while the taxi is deadheading: 

(a) A free taxi, while engaged in idle deadheading to Point 1, cannot be recalled to 

Point 2 to pick up a passenger. 

(b) A free taxi engaged in idle deadheading to Point 1 may be recalled to pick up a 
passenger. This requires that communication be available between the points of 
demand and the taxi. Since deadheading is required either before or after every 
call, the total running time per unit time in this case is F = 2v7 = 2x, 

To find the total waiting time under Policy A, we note that a call arising at Point 1 has 

to wait if the taxi is on its way for a previous call, in which case the average waiting time is 7; 
the total passenger waiting time at Point 1 is then vyFr. At point 2, the waiting time is 7 if 
the taxi is free and stationed at Point 1 (a fraction 1-F of the time); as before, the mean wait- 
ing time is 7 if the taxi is running (a fraction F of the time). The total passenger waiting 
time is, therefore, 


(6a) TR =U, FT + vo [((1-F)7 + Ft] = qx + 2px? ‘ 


Under Policy B (communications available), the probability is v/ V that the taxi is 
empty when traveling from Point 2 to Point 1. In this case, it may be recalled if a demand 
arises at Point 2; the average passenger waiting time is thus 7/2. The probability is Vo/ Vv 
that a taxi going from Point 2 to Point 1 is carrying a passenger, in which case the waiting 
time is (3/2)7. If the taxi is underway from Point 1 to Point 2 (probability F/2), the average 
waiting time for a passenger at Point 2 is 7/2. Hence, the total passenger waiting time is 


b F/"47.. "231 FT qa\.2 
6b Ta =, FT +v 1-F =(—< , A & xt @ 2p/1- y 
(6b) B 1 +a yr F(z, Bnet qx + P( 2) 
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Having communications in the taxi has shortened the total average passenger waiting 
time by pqx. In what follows, we assume that Policy B is in effect: the taxis may be 
recalled. 

Comparison. We now compare the total waiting times incurred under the two different 
idle-time strategies. This is done most easily by considering the ratio Tp/ fy A? where T A is 
the total passenger waiting time resulting under Policy A (Eq. (5)), and Tp is the total pas- 
senger waiting time under Policy B (Eq. (6b)). Instead of the parameters p and q, it is more 
convenient to use the ratio of the source demands r = (v,/ Vo) = p/q. The expression for the 
ratio of waiting times becomes 


(7) 1+r+r(2r + 1)x , 


Ta ter? 4—*2®_ (14 r+ r%)x 
(14+r)? 





Equation (7) is plotted in Figure 2 for different values of the ratio of source demands r, asa 
function of the utilization level x. From this plot, the following conclusions may be drawn. 
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Figure 2 - Ratio oftotal service time— scheduling 
and communication nonscheduling 


(a) As long as there is a real difference in source demand (i.e., r is 2 or larger) so 
that an idle-time policy which always returns the taxi to Point 1 actually represents 
an attempt to simulate the distribution of the demand, Policy B will provide a con- 

siderable improvement in passenger waiting times over Policy A, for low levels of 

utilization (x << 0.5). 
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(b) The advantage of Policy B, as compared with Policy A, decreases as the degree of 
taxi utilization increases to 50 percent. In other words, when the level of demand 
is so high that the taxi is kept busy just meeting this demand, it becomes less 
worthwhile to attempt redistribution or idle deadheading. 

(c) If the utilization is over 50 percent, it is disadvantageous to attempt redistribution. 


EXTENDED REGION OF DEMAND 
Dispatching Strategies 
We now consider a fleet of taxis which serves an extended geographical region with 
many sources of possible demand. We assume that the distribution of demand is uniform 
throughout the region. The most obvious dispatching strategy can be stated as follows: 
(a) If free taxis are available, the dispatcher sends the free taxi which is nearest to 
the call. 
(b) If no free taxis are available, he assigns the taxis as they come free to the waiting 
passengers in the order in which they placed their calls. 
This strategy, although very simple, is extremely inefficient. Following it, a taxi may have to 
travel the length of the area to reach its call, while in the meantime another taxi.is likely to 
come free much closer to the source of demand. This strategy may be modified by considering 
the two rules above as being applicable, for a given call, only to the taxis already free or 
coming free within a radius R of the call. This dispatching radius R may then be chosen so 
as to minimize the average waiting time for each call, by balancing the hazard of making R 
too small so that insufficient taxis will come free within R, against making R too large so that 
the caller may wait a long time while the taxi has to deadhead to him. The optimum dispatching 
radius and associated average passenger waiting time are calculated in (a) Dispatching Radius, 
which follows. 





Even with the inclusion of a "dispatching radius," the above strategy is not the most 
efficient one possible. A better strategy is to assign the taxi which will be available at the 
place of call first, independently of whether the selected taxi is occupied or free at the time 
when the call comes in. This strategy assumes excellent radio communication between the 
taxis and the dispatcher. Furthermore, a considerable degree of alertness is required on the 
part of the dispatcher, who must keep in constant touch with his taxis and must estimate for 
each occupied taxi the time at which the taxi will reach its destination and become available 
for further dispatching. This need for close communication was not present under the "dis- 
patching radius" strategy, when the dispatcher had only to keep track of free taxis. The addi- 
tional sophistication in dispatching results in much more efficient operations. In (b) The 
Optimum Dispatching Strategy, which follows, we derive the expressions for the distribution 
of passenger waiting times and number of taxis required to provide a given level of service 
under some general restrictions of exponential holding times and median waiting times small 
compared to the trip times, 

(a) Dispatching Radius 

Let A = area of region 





N = number of sources of demand 
n = number of taxis in the fleet 

7 = average trip time 

f = frequency of calls per source 
Vv = average taxi speed. 
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A call will either be served immediately, in which case the service time is the time it takes 
the taxi to travel to the caller—or the caller may have to wait for a taxi to become free, in 
which case the service time is lengthened by the time it takes to release a taxi. If there isa 
free taxi within R, this taxi will on the average be a distance r from the caller, where r is 
given by: 





(8) 


| 
iT] 
— 
yy 
w 
a) 
5 
w 
5 


«2 
= aR. 


The average time per assignment for a taxi is the average time of getting to the caller, 
T= 2R/3v, plus the average time of the actual trip, 7. Since the average number of trips 
originating per unit time is Nf, the average number of busy taxis is 


(9) n, = Nf (7 +28). 


Since there are (n - no) free taxis, the probability that there are no taxis in R is 


n-n 
7 R2 ' 
P, (R) = ( 7 re ’ 


i.e., the probability that all free taxis are outside R. The probability that there is a free taxi 
within R is 


9 oO 
i i oe “— « 2. 
(10) P, (R) = 1 Pa = 1 a . 
In the steady state, the average number of taxis released per unit of time must equal the 
average number of taxis picked up, Nf. The number of taxis released within R per unit time, 
therefore, is Nf(z R2/A). Hence, the average interval from a time of no free taxi in R till the 
time that the first taxi comes free is 


2 


(11) T, = A/27R°NE. 


The factor 1/2 expresses the fact that, on the average, a call occurs midway between succes- 
sive releases. The average passenger waiting time may now be written as 


n-n 


2 Oo 
(12) T = Poult, + 7) + Py, —— + 
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2 


If 7R° << A, Eq. (12) may be approximated by 





2 


2 
(13) r= —A_|1- (a) 138. 
27 R2NE A 3 v 


Substitution of Eq. (9) in Eq. (13) yields 


(14) 9 6 ente #5 th 
on R2NE Vv 2f N 


To find the value Rin which minimizes T, we differentiate Eq. (14) with respect to R, and 
equate the result to zero to obtain 


3 

yj 
15 R = ——-, 
ated m aNf 


The corresponding average waiting time 7. is found by substituting Eq. (15) in Eq. (14), 


yielding 
3 
/ & ‘ee 4). 
amNfv2 2f N 


Equation (16) shows the obvious result that we can make Ln smaller by increasing the number 
of taxis n. 


(16) <a 





bo | co 


For a numerical example, suppose that taxis with a cruising speed of 20 mph are to 
serve an area of 50 square miles in which an average of 40 calls arise per hour. The optimum 
dispatching radius is then 2 miles. If the average trip takes 0.2 hours and there are four taxis 
serving the region, then the average passenger waiting time is 0.2 hours. 

(b) The Optimum Dispatching Strategy 

In the following, we assume a call is always served with that taxi which can be 
snade available at the place of call first. If an occupied taxi is assigned under this policy, the 
total service time will consist of two parts: 

(a) the time required to complete the present trip; 

(b) the time required to reach the caller, starting from the point of release. 





Let 
n = probability that a given taxi is free = average fraction of free taxis. 
€ = probability that a given taxi is occupied = 1 - 7. 
The finite area of the region limits the part of the service time due to the travel of the taxi 
from its point of release to the location of the caller. Suppose that this limit is T O° Further- 
more, let an occupied taxi have a probability ~ (x) dx of coming free in a time interval 
(x, x + dx) from now. The probability density of service times Y(t) may be expressed in 
terms of the probability density of the continuing trip length, as follows: 
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(17) 


The terms appearing in Eq. (17) are contributed by: 
(a) The probability that the taxi is free, and that it is a time t away from the call. 
The probability of being a distance r away from the call is 27rdr/A. Since 
r = vt, the probability of being a time t away is Qnv2 tdt/A. It is clear that t 
has to be smaller than To: 
(b) The probability that the taxi is occupied, in which case the probability that it will 
come free in a time x is ¢(x)dx, and the probability of having an additional time 
(t- x) to goto the call is 27 v2 (t - x) dt/A. 

Since only the total service time is of interest, it is necessary to integrate over the 
possible continued time of the present trip. This integral runs from t- T> to t, unless t is 
smaller than To in which case the interval is from 0 to t, since x cnante be negative. If 
V(t) is to be iit normalized, i.e., Sg y(t) dt = 1, t must be taken as the mean time- 
radius of A, so that 


ber 


Let X(t) be the probability density that the taxi with the least service time is time t 
away. This is then also the waiting time distribution for a caller. X(t) may be expressed in 
— terms of W(t): 


xis 
t n-1 
(18) X(t) = nv¥ (t) - { vera ‘ 
(0) 
1e 
This is because 
t 
{ y(t") dt" 


(0) 
is the probability that the service time is less than t for a given taxi, so that 
r- t ° 
: y(t") dt 
0 


is the probability that it is longer than t. The factor n results from the fact that it is imma- 
terial which of the n taxis is the closest. For n >> 1, 
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t 
(19) X(t) = nw (t) exp a | W(t") dt"| . 


If, for example, ¢ (x) = (1/7) e7X/ " (exponential holding times) then from Eq. (17) we obtain: 


2 2 
_ 21Vv ét 
(20) y(t) = — (ne +4 2) fort << T. 


Equation (20) is valid if the passenger waiting time is generally less than the average trip 
time. In this case Eq. (19) becomes: 


Pit et fe¥ 
(21) X(t) = a [nt + ; () exp fag 2 


where a, = Qnv" 72/A, 


Taxi Fleet Size and Service Time 

The size of the taxi fleet has to be chosen sufficiently large to provide the desired level 
of passenger service. The'criterion we take for the desired level of passenger service is that 
all but a fraction € of the calls have to wait less than a time t; in other words, the probability 
that a passenger has to wait longer than t is €. This includes as a particular case, obtained by 
setting € = 1/2, the criterion that t be the median waiting time. 

From Eq, (21) the probability of a waiting time longer than t is given by 


2 3 
(22) P (waiting time >t) = exp {om F n (£) + : (£) }} 


where 7 is the mean trip time. 





For mathematical convenience we measure t in units of 7, calling a = t/7, so that our 
service criterion becomes 


(23) P (waiting time = aT) =e. 


Inserting (23) into (22) and rearranging, it is found that: 





(24) n= nn(1- 3) Sine, a <<1, 
a 
age 
Equation (24) constitutes one relation between the number of free taxis (n7) and the total 
number of taxis n, 

Another relationship between n and (nn) may be obtained from the requirement that 
the total occupied time of the taxis must equal the number of trips made, times the average 
time spent per trip. The average taxi trip time consists of two contributions: the average 
deadhead time 
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t = t X(t) dt , 


plus the average time of the trip the passenger desires to make, 7. If N potential passengers 
each make an average of f calls per hour, this requirement results in the equation, 


(25) né = n(l - yn) = Nf(t +7). 


If t is computed in the approximation that most calls are served with free taxis, 


a.n7t 2 

X(t) = = exp |-a njt 
2 ° 3 

T aT 








and 


o aa 
2v nn 


| 


(26) 


Substituting Eq. (26) into Eq. (25) and solving for n, we obtain 


(27) n= nn+nt Asner. 
Vv 


n?) 


Equations (24) and (27) are two simultaneous equations for the total number of taxis n and the 
average number of free taxis n7, resulting from the requirement that € of the calls have to 
wait a7 or longer. Equating (24) and (27) yields one equation for n: 





3 ao? 6v 


(28) (nn)3/2 + (My | (nn)!/2 , Nia ae 

The total number of taxis required to meet the service criterion can be obtained by using the 
solution of Eq. (28) in Eq. (24). The distribution of service times can be obtained from Eq. (21). 
Letting € = 1/2, these equations serve to set the median waiting time at a7. Since the approx- 
imations are valid only for a << 1, these equations are valid only in the case of median wait- 
ing times which are small compared to the average trip time. 


MONTE CARLO TAXI SIMULATION 
Where there are several sources and taxis and demand is randomly distributed, it is in 
general too difficult to develop optimum scheduling rules and to calculate the passenger waiting 
times and taxi deadheading by analytic means. In order to test the conclusions developed by a 
consideration of more idealized models (presented earlier in this article), a Monte Carlo 
simulation was carried out for a closed-system model involving four source and destination 
sites and eight taxis. The sites are geographically distributed on the corners of a square 
(Figure 3). The calls for taxis come randomly at an average of 1, 2, 2, 3 per hour for Sites A, 
B, C, D, respectively. The destinations are completely random. The trip times are .35 hours 
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between sites A-D and B-C, and .25 hours between all other 
sites. A taxi "stand" is located at the center of the square, 
distant .175 hours travel time from any source site. 











ra 2 Three dispatching policies for utilizing the free taxi 
© 6 time were examined in the simulation: 
S “ (a) Leave the taxis at the destination point until a 
demand arises. 
(b) Send the taxis to the central stand as soon as they 
0.25HRS become free. 
c D 


Figure 3 - Distribution 
and 


of s 
taxis 


demand first. It was assumed that communications between scheduler, sites, and taxis are 
excellent, and in Case 2, that a taxi could be redirected with no loss of time. 

We assume that the number of calls per hour is distributed in a Poisson distribution. 
It follows that the distribution of time intervals t between calls is exponential: 


where m = 1, 2, 2, 3 is the mean number of calls per hour for sites A, B, C, D, respectively. 
This exponential relation was used to set up a random demand schedule for each site. The 
destinations were chosen randomly with equal probabilities. 

The results of carrying through the simulation for the equivalent of a 9-hour period are 
summarized in Table 1. The distribution of passenger waiting times is presented in Figure 4 
for each of the three cases. In Case 1, the taxi deadheading was always in answer to a call 
from another site. The deadheading in Case 2 was mostly to and from the stand; in a few 
cases, the taxi was required to proceed to another site immediately after discharging a pas- 
senger. In Case 3, the deadheading turned out to be almost entirely for redistribution purposes 
to meet possible demands. 


ites, demand, 


(c) Redistribute the free taxis among the four sites in 
a ratio as close to 1:2:2:3 as possible. 

The taxi stand is used in one case only, Case 2. For all cases, 
whenever it was necessary for a passenger to wait for a taxi, 


the taxi was assigned which would arrive (free) at the point of 


P(t) = — exp - (t/m), 











TABLE 1 
Results of Simulation, Fixed Number of Taxis 
Total No. Passenger Passenger Taxi Taxi 
Case of Travel Time | Waiting Time | Deadhead Time | Running Time 

Calls (hr) (hr) (hr) (%) 
1 86 23.70 6.21 5.90 40.5 

2 86 23.70 11.84 15.83 55 

86 23.70 1.33 7.20 43 
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Figure 4 - Passenger waiting times--Cases 1 through 5 


In Case 1 (no redistribution) the total passenger waiting time was 6.21 hours; the taxis 
were forced to deadhead for 5.50 hours. The result of using a stand in Case 2, as compared 
with Case 1, was to double the passenger waiting time and triple the taxi deadheading time. All 
passengers were forced to wait under this policy; the peak waiting period was shifted, however, 
toa smaller waiting time. Thus, although the taxi running time is increased in Case 2 to 55 
percent as compared to 40.5 percent in Case 1, very little advantage was accrued. 

In Case 3 (redistribution), the total passenger waiting time was reduced from 6.21 hours 
to 1.33 hours, with an increase in deadheading of 31 percent over Case 1. With this moderate 
increase in deadheading for redistribution, the mean waiting time per passenger was reduced 
80 percent, using the same number of taxis. The advantage of good redistribution policies are 
thus markedly exhibited. 


Variable Number of Taxis 

It is interesting to extend the Monte Carlo analysis to examine the effect of the number 
of taxis. The same geographic model and demand schedule were used for the case where only 
six taxis (Case 4) or only four taxis (Case 5) were available. Wherever possible, redistribu- 
tion was made. The initial conditions for the simulation are listed in Table 2. The initial con- 
ditions for Case 3 (eight taxis available) are included to complete the comparison. 

The results of the simulation are summarized in Table 3, and the distributions of pas- 
senger waiting times are plotted in Figure 5. In Case 4, 80 percent of the deadheading was for 
redistribution; in Case 5, the taxis were utilized so fully (74 percent of the time) that there was 
little opportunity for redistribution. Both the total taxi deadheading and total passenger travel 
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TABLE 2 
Initial Conditions, Variable Number of Taxis 

















No. of Taxis 
Site 
Case 3 | Case 4 | Case 5 
A 1 1 1 
B 2 1 1 
Cc 2 2 1 
D 3 2 1 
Total Taxis Available: 8 6 4 




















TABLE 3 
Results of Simulation, Variable Number of Taxis 











No. of Total Passenger | Total Passenger Total Taxi Taxi 
Case Taxis No. of | Travel Time Waiting Time Deadhead Time | Running Time 
Available | Calls (hr) (hr) (hr) (%) 
3 8 86 23.70 1.33 7.20 43 
4 6 86 23.70 5.01 7.75 58 
5 4 86 22.35 11.98 4.10 74 


























time were reduced (due to the combining of trips) in Case 5. As expected, the mean waiting 
time per passenger is increased from 0.9 minutes in Case 3, to 3.5 minutes and 8.4 minutes 
for Cases 4 and 5, respectively. In addition, the maximum service time under Case 5 was 
increased to 0.6 hours, as compared to 0.3 hours under Case 3. 

The simulation confirms the results obtained from the theoretical models: Redistribu- 
tion to match expected demand decreases the passenger waiting times, but this effect becomes 
less important as the level of taxi utilization is increased. 
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ON NETWORK FLOW FUNCTIONS*t 
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This article covers an investigation of the capacity of a network as a 
function of the capacities of its individual arcs. The case of two vari- 
ables was investigated in detail, and it was found that each pair of arcs 
either consistently help each other or consistently hinder each other. 
Interaction types were computed for a number of special cases, such 
as arcs in parallel or arcs in series. 











7] INTRODUCTION 


This article is concerned with the structure of the flow through a capacitated network. 
Formulas are derived for the maximum flow through the network as a function of the individual 
arc capacities. The two-variable case is presented in detail, and it is shown that certain pairs 
of arcs stand in a complementary relationship to one another, whereas other pairs are mutually 
interfering. Moreover, the interaction type of a given pair of arcs is the same for all values of 
their capacities. In a number of special cases the interaction type depends only on the topological 
structure of the network, and is independent of the numerical capacity values. 


THE MAXIMUM FLOW FUNCTIONS 

Let I be a directed graph with stated, nonnegative capacities Cy, Co, +++» Cy ON its 
arcS @), A, «++, A. Certain nodes of I, which must have only out-pointing arcs, are 
designated as sources; certain others, with only in-pointing arcs, are designated as sinks; the 
a remaining nodes might be called intermediate. Let " denote the network so defined. A flow 
S pattern on {i may be described by a nonnegative vector f = (f,, eke fy) < (c,, eieteus Cc.) such 
that at all nodes other than sources or sinks the sum of the fi on in-arcs is equal to the sum on 
out-arcs. The flow through fis then defined as the sum of the f; over all out-arcs at sources 
(or, equivalently, over all in-arcs at sinks). Let 3(f) denote this sum, and let F denote its 
maximum over all flow patterns f; this is sometimes called the “capacity” of . We propose 
to study F as a function of the cj. 

It will be convenient occasionally to allow infinite arc capacities, and hence to allow F 
to be infinite. However, we assume that the originally given c, are finite. 
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Let ni (x) denote the network obtained by replacing the capacity ci of @ i by the non- 
negative variable x, holding everything else fixed. Let F; (x) be the capacity of ny (x). Then 


we have 
F; (0) +x 


(Al) F, (x) = min 
F; («) 


« 
In other words, for x below a certain critical value, namely, x = Fi (co) - Fi (0), we have 
i. 
dFj/dx = 1, whereas for x > x we have d F,/dx = 0 (see Figure 1). 











Figure 1 


Similarly, if ; and a, are distinct arcs of T', we may define F;; (x, y) to be the 
maximum flow through Nii (x, y), the network obtained by replacing ci by x and cj by y. 
Then we have 


F,, (0, 0)+x+y 

F..(0, 0) + x 

F.. (x, y) = min 4 
1) Fi (0, 0) + y 


F ij (~, 00) 
which follows from (A1) by virtue of 


Fij (0, y) +x min [Fi (0, 0) + x+y, Fi; (0, 0) + x] 
Fi, (x, y) = min 
Fij (co, y) min [Fi (co, 0) + y, Fij (2, co) | 


By induction, similar expressions (A3), (A4), ... , can be derived for the functions 

Fig % y; 2), Fike ee, sce, « OB we carry this procedure through to the end, making 
all arc capacities variable, we see that the Fio n& ..) term of the minimand of (An) is 
always either 0 or o, and the following expression results: ° 


Fi... n (Byers > Xp) = min » Xi» 
S ieS 
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where the "min" runs over all subsets S of {1, 2,...,n} whose complements S' contain no 
chain from a source to a sink. The reader will recognize this as a form of the "min cut - max 
flow’ theorem.! 


THE TWO-VARIABLE CASE 

We shall now investigate the possibilities for the function Fi, (x, y) in more detail. 
Referring to (A2), we see that any such function will partition the positive quadrant of the x, 
y-plane into at most four open, convex regions in which it is linear, plus certain boundary lines 
and vertices. We label these regions according to the values of the partial derivatives: 


OF ,;/0x aF,,/ey 





1 1 





1 0 





0 1 





0 0 














Continuity of Fij requires that the common boundaries of each pair of regions be oriented in a 
definite way, as shown in Figure 2. 


R 
he * 00 
RiilFor = RA AR, Bao Foo a 


(b) (c) (d) (e) (f) 


Figure 2 


Thus Roo» for example, is unbounded to the top and to the right, etc. In fact, there are only 
three possible "complete" configurations, as shown in Figure 3. 

The numerous "incomplete" cases that occur when one or more of the regions are 
missing prove to be subconfigurations of the ones illustrated. We note that in all cases there 
is never more than one diagonal boundary segment, since an Rii - Roo contact precludes an 
Rio - Roi contact. 

We are now in a position to make a general statement about the behavior of Fij (x, y). 
First, we observe that because of the piecewise linear nature of Fij the second derivatives 
are not particularly informative; they are always either zero or infinite. However, we can set 
up the finite difference quotient: 





Isee L. R. Ford and D. R. Fulkerson, ' Maximal Flow Through a Network,'' Canadian Journal of 
Mathematics, 8 (1956), pp. 399-404. We refrain from claiming a new proof of this celebrated 
theorem, since the argument is based on statement (Al). Although an immediate consequence 
of the theorem, (Al) is not particularly easy to establish from scratch. 
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(b) 
Figure 3 


2 
4 A Fi, : Fi (x+h, y+k) - Fi (x+h, y) - Fij (x, y+k) + Fi (x, y) 


Si Ax Ay 7 hk : 





in place of a2 F,,/0x dy. This is of course a function of h and k as well as x and y, and is 
well-defined only if h = - x, k =- y, and hk ~ 0. 


THEOREM 1: Within the class of networks ni (x, y), x = 0, y = 0, the quotient Qi; is 
always = 0, or is always = 0. 


REMARK: Roughly speaking, a positive Qi means that the arcs a; and @. complement 
(reinforce) each other; a negative Qi; means that they compete (interfere) with each other. 
Typical examples are arcs connected "in series," and arcs connected "in parallel," respec- 


tively. (See the next section.) 
PROOF OF THE THEOREM: Without loss of generality, we may assume h >0,k> 0. 


Consider the rectangle with vertices (x, y), (x +h, y), (x +h, y +k), (x, y +k). It cannot enclose 
more than one diagonal piece from the boundary configuration as defined above. If it encloses 
none, then we immediately have Q; = 0. If the diagonal piece enclosed has positive slope 
(Figure 2(d)), then Qi; is strictly positive; in fact, it is exactly equal to the length of the inter- 
cepted diagonal segment divided by V2 hk. Conversely, a negative slope (Figure 2(c)) yields a 
strictly negative Qi. These two cases being mutually exclusive, the theorem follows. 

Letting the rectangle of the preceding proof expand to infinity, we obtain the following 
corollary, which relates the signs of the Qi;'8 to the signs of certain constants Li}, given by: 


Li; = Fi (co, co) - Fi, (co, 0) - Fij (0, 0) + Fi, (0, 0) , 


which, in a given network, are generally easier to evaluate than the Qi;'8- 





COROLLARY: (a) If Lij > 0 (including Li = +0), then Q: is sometimes > 0, and 
never <0. (b) If L,, = 0, or if Li is indeterminate (© - ©), then Qi; is always = 0. 
(c) If Li; <0, then Qj is sometimes < 0, and never > 0. 
















ON NETWORK FLOW FUNCTIONS 



































SOME SPECIAL CASES 
Theorem 1 enables us to classify the pairs of arcs of a network Vi as either comple- 
mentary (Q;, sometimes > 0), interfering (Q;, sometimes < 0), or neutral (Q;, £0). In 
general, the interaction type of two arcs will depend on the capacities of the other arcs of fl, 
as well as on the relative positions of the two arcs in the network topology.* In this section we 
exhibit several cases in which topological considerations alone suffice to determine the inter- 
action type, independently of the numerical capacity values. The results may be summarized: 
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Interfering Complementary 





Arcs in parallel Arcs in series 
Confluent arcs 
Divergent arcs Source - sink 
Source - source 
Sink - sink 














First, we shall state and prove two simple lemmas. A cut can be defined as a collection 
of arcs that intersects every (directed) chain from source to sink. The "capacity" of a cut is 
the sum of the capacities of its arcs; a cut of minimum capacity will be called, for brevity, a 
S min cut. The set of nodes that can be reached from the sources without using the arcs of a 

given cut will be called the source set of that cut. 
ent 


LEMMA 1: If A and B are the source sets of two min cuts, then the cut consisting of 
all arcs from nodes in A  B to nodes not in Af B is also a min cut. 


PROOF: Let F be the maximum flow through ii. A flow pattern f achieving that 


0. maximum must saturate the arcs of all min cuts, and cannot include any "back flow" into their 

ose source sets. Hence the flow along the arcs leading from A / B to its complement will total F, 

5 and will saturate these arcs, These arcs therefore constitute a cut of capacity F — the mini- 
mum value. 

ph. 

a LEMMA 2: If the capacity of an arc as is at the "critical" value 


 * F,(«) - F,(0) 


(see Figure 1), then that arc (a) belongs to at least one min cut, and (b) fails to belong to at 
least one min cut. 


PROOF: Decreasing the capacity of ar decreases the maximum flow through /; hence 
(a) follows. Increasing the capacity of @ i does not increase the maximum flow; hence (b) 
follows. 








The example in the last section of this article is a case in point. 
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THEOREM 2a (ARCS IN SERIES): If the terminal node of a, is the initial node of a, 
then Qi; is always > 0. 


PROOF: Let x and y be positive, and consider 
a min cut in Nis (x, y). If the common node, N, is the 
source set, then a@ i cannot belong to the cut; if N is not 
in the source set, then @. cannot belong to the cut. 
Since Rij is the region in which both arcs are in every 
min cut, we see that Rij is empty. Hence Qi; =U, 


THEOREM 2b (ARCS IN PARALLEL): If as and *, have the same initial node and 
the same terminal node, then Qi; is always = 0. 


PROOF: This is a special case of divergent or 
confluent arcs, as follows. 


THEOREM 2c (DIVERGENT ARCS): If a; and a, have the same initial node, then 
is always = 0. 


Qi; 


PROOF: If region Rp, is empty, then Fij (x, y) 

is unbounded, implying that at least one of a p a; isa 

direct source-to-sink connection, and giving us Qi = 0. 

Therefore we assume Roo nonempty. Let (x*, y*) be a 

minimal point in its closure.? Since the theorem is 

immediate if either Roi or Rio is empty, we can 

require that both x* and y* be strictly positive. Hence, 
in the network Maj (x*, y), both a; and a. are critical, in the sense of Lemma 2, and both 
belong to min cuts. The source sets of these cuts (which may or may not be identical) must 
both contain the common initial node of a, a. Moreover, neither terminal node can belong to 
both source sets. Therefore we can apply Lemma 1 to the two cuts to obtain a single min cut 
to which both arcs belong; the capacity of this cut is of course Fij (x*, y). In the variable 
network Nj (x, y), its capacity is then: 


Fij 


x * * + 
(x,y)+x-x +y-y, 


which gives us an upper bound for the max flow Fij (x, y). In particular, we have 


x * * * 
Fj; (0, 0) < Fj x +2 ©¢9 





3 That is, a point such that x < x* 


and y < y* never hold simultaneously for any other point 
(x, y) in the closure of Rgo. 
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By (A2), equality holds here, and we see that the point (x", y) must be in Ri or its bound- 
ary. Thus, Roo and Ri touch, and we have Qj = 0, as required. 


THEOREM 2d (CONFLUENT ARCS): If as and a@ j have the same terminal node, then 
Qj is always = 0. 


PROOF: Obtained from the preceding proof by reversing 
the direction of the network. 


THEOREM 2e (SOURCE, SINK): If the initial node of a, is a source and the terminal 
node of *, is a sink, then Qj is always = 0. 


PROOF: If Roo } is empty, the result is immediate. If Roo 
is not empty, then let (x* >y *) be a minimal point of its closure. 
Excluding the trivial, neutral case where no diagonal segments 
exist anywhere in the boundary configuration, we may demand that 
both x andy _ be strictly positive, and hence critical. Lemma 2 
then gives us a min cut to which or belongs and also a min cut 
(possibly the same one) to which a; does not belong. Applying 
Lemma 1, we obtain a single cut whee source set contains the initial node of a, (since it isa 
source) but not the initial node of @. (since it is connected to a sink), nor the wunied nodes of 
either. arc. In short, the cut contains a; but not aie Its capacity in general is equal to 
Fi s(x" wy )+x-x- . Thus, for all y, 


* * * 
Fij (x, y) Fj, »y ) +xX- X 
Symmetrically, for all x, 
* * * 
Fij (x »y ) +7 > = 


In particular, for all h, 


* * * * 
F jj +h,y - h) aa Fj, (x »y ) - |h|. 


It follows that a diagonal segment in the minimal set of Roo (Figure 3(c)) is impossible. 
Therefore Qi; = @, 


THEOREM 2f (SOURCE, SOURCE): If the initial nodes of @; and *, are both sources, 
then Qj is always = 0. 


PROOF: Reduces to 2c. 
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THEOREM 2g (SINK, SINK): If the terminals of a i and a j are both sinks, then Qi; 
is always < 0. 





PROOF: Reduces to 2d. 


AN EXAMPLE 

One might be tempted to generalize the idea of "parallel" arcs beyond the rather trivial 
case enunciated in Theorem 2b. Thus, in the network Ny (x, y) shown in Figure 4, @ 1 and a, 
are "parallel" in the sense that no directed chain contains them both. However, their inter- 
action is of the "serial" rather than the "parallel" type, since we have Qi =+1 for 
x=y=h=ke=1, 











5 
>—® 
x=ye=l 
x=l,y=2?-F=4 
x=2 = 1 
5 » y 
ees K=y=2 -F=5 
Figure 4 


Curiously, removing the arc a, (or reducing its capacity to zero) has the effect of 
inverting the interaction-type of A1, %, since then x = y= 2, h=k=1 yields Qi =-1, 
This latter variation suggests a definition of "parallel chain," which would lead to a minor 
extension of Theorem 2b. But it seems evident that substantially more powerful topological 
criteria than those of Theorems 2a-2g are not to be expected. 
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PRINCIPLES OF LOGISTICS - A PROVISIONAL DEFINITION* 


R. M. Shoemaker, Lt Col, USAF 


School of Logistics, Air University 
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EDITOR'S NOTE: This article is unusual in that its form and content 
appear to be aimed at Military ''Commanders" rather than Logistics 
"Researchers' who in general are interested in specific techniques 
rather than foundation principles. However, the article covers a most 
important subject, and it is presented as an ''expository" article by a 
"Logistics Thinker" to stimulate greater interest in the field of 
Logistics. 











INTRODUCTION 

In the constant endeavor to develop techniques and systems through which the forces of 
military organization can be supported with true effectiveness, it is inconceivable that a body 
of principles to guide the endeavor may not exist. Yet, such appears to be the case with 
respect to the artful science of Logistics. One may argue (a) that principles do in fact exist— 
however obscure—in the field of logistics and (b) that certain of these principles have guided 
the development of military logistics systems into their present structures. This follows as a 
logical consequence of the fundamental purposes underlying the very existence of logistics: to 
assure that the quality and quantity of resources needed by military forces are in the hands of 
the forces in the time and condition needed to successfully fulfill the operational mission. The 








foregoing has been generally recognized and has undoubtedly served directly and indirectly as 
the basic guide in developing logistics systems and structures; it is, therefore, the basis for 
all logistics principles. 

The fact is not that principles do not exist, or are not followed in the field of logistics; 
it is, rather, that recognition of such principles is too much taken for granted. The principles 
are taken for granted to the extent of not being written down, much less generally understood. 
Therein lies the danger for, while adherence to principles is no positive assurance of success- 
ful venture, actions taken within the framework of a set of principles are more often right than 
wrong. But, the key to truly effective practice of logistics is not within any individual principle 
or collection of principles. The key is, rather, in the understanding of the relationships 
between the principles themselves and between the principles and the purpose(s) for which 
the very endeavor (logistics) exists. The successful use of the principles lies in knowing which 
principles apply with what order of emphasis to produce the most effective results in a given 
Situation. This is an aspect of application which follows comprehension of the principles 
themselves. This comprehension is the intended purpose of this article. 





*Manuscript received April 1960. 
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TOWARD A DEFINITION OF LOGISTICS PRINCIPLES 


Principles governing any field of endeavor are defined only after lengthy study and 
analysis of a whole group of actions taking place in that endeavor. The principles of war—as 
originally identified by Jomini and Clausewitz, for example—were deduced only after thorough 
study of a whole array of campaigns and battles back through the ages. 


There is one significant fact which may hasten the definition of logistics principles. 
The relationship between strategy and tactics (combat) and logistics (the support of forces 
both in combat and preparation of combat) is so completely intertwined that the principles 
which guide one must be discernible in principles which guide the other. "There are five 
things from which a soldier should never be separated: his musket, his cartridges, his 
knapsack, rations for at least 4 days, and his pioneer tool. Let the knapsack be reduced to the 
smallest size possible, if this is deemed necessary, but he should always have it with him."! 

There is, then, this means for deducing "Principles of Logistics": Strategy and tactics 
(combat) comprise the art and science of applying resources to accomplish military—and 
national—objectives; logistics comprises the art and science of preparing resources for 
application through strategy and tactics. The relationship of strategy-tactics-logistics is, 
therefore, so directly meshed as to produce quite similar principles for the conduct of all three 
efforts. "Principles of War" are widely recognized and thoroughly tested. These should pro- 
vide a justifiable basis from which to deduce "Principles of Logistics.'"' Principles of logistics 
must complement, supplement, and, in some cases, be nearly identical with principles of war. 

On this belief, the Principles of Logistics—admittedly provisional until further studied— 
are defined as: (a) The Principle of the Objective; (b) The Principle of Interdependency with 
Strategy and Tactics; (c) The Principle of Coordination; (d) The Principle of Support Integra- 
tion; (e) The Principle of Evolution; (f) The Principle of Dispersion; (g) The Principle of 
Flexibility; (h) The Principle of Control; (i) The Principle of Economy of Resources; and 
(j) The Principle of Personnel Participation. These principles are discussed individually in 
the remaining part of this article. 


PRINCIPLE OF THE OBJECTIVE 


Every logistics endeavor must be guided by a clearly stated 
objective. The objective, to the greatest extent possible, must be so 
specified as to permit continuous measurement of the degree of 
accomplishment of the endeavor toward the objective. The logistics 
objective is valid only insofar as it supports the combat force objective; 
therefore, the propriety of the logistics objective must constantly be 
reappraised in the light of the intentions of the combat force. 


Logistics, as a component of the total military effort, has but one fundamental reason 
for being: to support the combat forces. To provide this support in the most effective manner, 
still assuring reasonable economy of resources, becomes the true objective of the logistics 
effort. This total objective of higher level logistics endeavor may, and must, be translated 
into component subobjectives for the individual command, the base, the unit, and the individual. 





1Napoleon, Maxim LIX. 
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The efforts of each must be credited as contributing to the whole. This can be the case only to 
the extent that each is given an objective, recognizes this objective, and works toward this 
objective. There can be no greater assurance of effective logistics than the existence of a 
firm sense of purpose in all who are responsible for logistics. This sense of purpose can only 
be defined in terms of the objective, that for which logistics exist. 

Logistics objectives are properly established in two mutually dependent categories: 
the ultimate objective(s) and the attainable objective(s). The attainable objective exists within 
the structure of the ultimate objective, but the ultimate objective is pursued as a long-range 
endeavor only through fulfillment of the intermediate, attainable objective. Ultimate objectives 
may or may not be attainable; they may exist essentially to give sense of direction and motiva- 
tion to the logistics effort or they may be truly attainable, but only over an extended period of 
time. Attainable objectives, conversely, recognize the constraints of resources and are 
accordingly developed as successive steps toward the ultimate objective. In a sense, ultimate 
objectives are of a strategic nature, attainable objectives of a tactical nature. 

The objective assumes particular significance in relation to the elements of integration, 
coordination, and control inherent in the pursuit of logistics. Logistics is characterized by 
diversification of organization, functions, and resources. Effective performance, therefore, 
requires a compensation for this diversification; the objective and its accompanying sense of 
purpose provide this compensating factor. The objective serves as a common guide to the 
various activities which perform logistics; most important, it serves as a cohesive agent 
which tends to pull this diversification back together. The objective, then, is of paramount 
importance to the logistics effort. Clear recognition of the objective and acceptance of the 
objective as the guiding influence in any logistics endeavor make possible the most positive 
assurance of effectiveness. 

There are three distinct divisions to the fundamental logistics objective of force 
support: to establish, to sustain, and to modernize. Of these, the third must be given priority 
during peacetime; the second during combat. 

Examples 

Certain examples may serve to illustrate the application of this logistics principle. 
Within the Air Force, the Air Force Logistics Command, formerly the Air Materiel Command, 
for many years has used the slogan “More Air Force Per Dollar.” Whether so recognized or not, 
this statement actually existed as an objective, an ultimate objective toward which all of the AFLC 
efforts were directed. Who was to say when this objective was attained? It is obviously a matter 
of conjecture, yet, the objective has served an extremely useful purpose as a “direction pointer,” 
a motivation, for the over-all efforts of AFLC in the materiel aspects of logistics. 

Numerous attainable objectives have been set by the AFLC within the framework of this 
ultimate objective. For certain selected aircraft, a goal of zero was set for the rate at which 
these aircraft existed in an "out-of-commission" status; this obviously contributed directly 
toward the "More Air Force" portion of the ultimate objective. Each of the units of AFLC 
responsible for ''in-commission" status of aircraft directed their individual efforts toward the 
zero objective and thus contributed to the ultimate command objective. A related attainable 
objective exists in the area of selective management of Air Force materiel. Briefly, the 
Air Force recognizes that its inventories can be grouped according to the costs of materiel 
involved and that different emphases of management are warranted by the different groups. 

In this fashion, greatest management attention can be devoted to the relatively few items which 
represent a major portion of the AF inventory dollar. Asa result of these different management 





162 R. M. SHOEMAKER 


emphases, it has been possible to establish a series of specific inventory stockage objectives 
(expressed in number of days of supply) throughout the world which represent a minimum 
level without jeopardizing customer support. These attainable objectives are on the "per 
dollar" side of the over-all command objective—that of directing individual unit efforts toward 
a specific goal. 


PRINCIPLE OF INTERDEPENDENCY WITH STRATEGY AND TACTICS 


Logistics must be interdependent with strategy and tactics, since 
logistics involves the preparation of resources for application through 
strategy and tactics by the combat forces. The dependency of strategy 
and tactics upon the performance of logistics and the dependency of 
logistics upon the planning of strategy and tactics require recognition 
and application at all levels and by all personnel. 


The strength of military forces and the means of combat exist in resources. There is 
a preparation of these resources which is logistics, and there is an application of these 
resources which is strategy and tactics. Through a combination of this preparation and 
application, military force is exerted to achieve objectives. Key to this exertion is effective 
combination, a situation which exists only to the extent that the interdependence of strategy and 
tactics and logistics is recognized and applied in the conduct of all three. 

Recognition of this interdependence is rooted in acceptance of the fact that the elements 
of military operation and administration, and the command context, in which both exist, are 
inherent equally in strategy, tactics, and logistics. That is to say, this principle underscores 
the existence of "logistics operations" as well as "tactical and strategic operations," "logistics 
command" as well as "tactical and strategic command." "Operations" and "command" cannot 
be construed to exclude logistics, since logistics is their very life blood. Just as there are 
strategists and tacticians, there are logisticians. Also there are tactical operations and 
logistical operations, and there are command of tactical and strategic efforts and command of 
logistical efforts. Artificial barriers which tend to ignore the mutuality of strategy, tactics, 
and logistics jeopardize the true relationship. 

The relationship is not difficult to visualize. The progression, always starting with the 
objective, runs: objective—strategy—tactics—combat techniques—forces—logistics. Logistics 
cannot be placed at the head of the progression, since it cannot rightly specify any of the other 
elements of the chain. Still, it must be an integral part of the progression, since without 
logistics, the progression is impossible. Objectives, strategies, tactics, and the like—which 
are developed independently of logistics capabilities—have the mark of failure. "For want of a 
nail - - -"' Logistics fits properly at the end of the progression, at the point where accomplish- 
ment of the objective is detailed in specific forces and associated tasks. These forces and 
tasks, translated into resource needs, are the instigators of logistics. The extent to which 
logistics fulfills the resource needs, as this extent passes upward through the progression, 
bears significantly on achievement of the objective. 

Logistics asks one fundamental question: What forces are to be supported, where, and 
under what conditions? Strategy and tactics answer. It is in the context of strategy and tactics 
that forces are determined. It is the forces, not the strategy or tactics employed, which 
logistics supports. 





















































PRINCIPLES OF LOGISTICS 


To fail to recognize the interdependency of strategy—tactics—logistics is to fail to 
grasp the fundamental principle of warfare: "To utilize strategic maneuver to bring the 
a greatest mass of one's forces in a coordinated effort upon a decisive point of the enemy's 
rd 4 forces.""~ Adherence to this principle is scarcely probable, lacking effective integration of 
i logistics, tactics, and strategy. It is, conversely, entirely possible when the three are meshed. 
Initial planning of a strategy and tactics is the proper place for application of this 
interdependency principle to begin. Planning which neglects it is suspect from the start. 
Examples 
That soldier who gave as his means of success—"'to get there first with the most"— 
give us an excellent though abbreviated illustration of this principle. Putting his ideas another 
way, his strategy called for early positioning of his forces which he endeavored to have in 
overwhelming numbers as compared to the opposition. His tactics involved, essentially, rapid 
troop movements. Both his strategy and his tactics depended upon the logistics of his forces: 
the means of movement, the daily provisioning of the force, and the resupply of arms and 
ammunition. His logistics were obviously dictated by his strategy and tactics of the day, and 
is vice versa. None of the elements could be satisfactorily successful without due regard for the 
others. 
In a more modern context application of this principle of interdependency is clearly 
evident in the "island-hopping" techniques of Gen. McArthur's forces in the Pacific during 
and World War I. It was quite evident that success of the basic strategy and specific tactics 
employed there was heavily dependent on the continuous effective performance of the logistics 


AW 


nts forces. These latter forces, to perform effectively, had to be advised, as an integral element 
of the campaign, of the forthcoming operations. 





Hics PRINCIPLE OF COORDINATION 
ot 
Coordination is the most essential ingredient of every logistics 
endeavor. It is the means whereby the logistics endeavor and the 
of forces to be supported are joined; it is the method by which individual 
segments of the total logistics effort are joined in a common endeavor. 
Support, lacking coordination, breeds cross purposes and loses objectivity. 


he 

cs Logistics, historically, has involved a manipulation of shortages in an effort to support 

er military forces. That this condition exists today is apparent to all. The logistics challenge is 
to do more with less. The ability to meet this challenge lies in coordination of efforts, both 

| between forces and logistics and between elements of logistics. 

fa The principle of coordination is closely related to another logistics principle, that of 

ish- support integration. Coordination involves the administrative-type process of assuring com- 
mon action among the various distinct elements of logistics. Integration involves a physical 
uniting of related elements so as to reduce the number of distinct elements which, while they 
exist as distinct elements, tend to complicate the total logistics process. As greater integra- 

tion is achieved, lesser coordination is required, and vice versa. The two principles are 
d accordingly inversely proportional in application. 
ics 





“Treatise on Grand Military Operations, Gen. Antoine H. Jomini. 
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Coordination is inherent in many of the other principles of logistics and, yet, is 
deserving of stature as a principle by itself. This is so because coordination is the very 
foundation of effectiveness and economy. It is, in a sense, both the means to and the result 
of efficiency. It is the key to flexibility. 

Coordination is not an act which should be looked upon as a crutch for poor organiza- 
tion or for ineffective control. It should, instead, stem naturally from the specific missions 
assigned various components of logistics. This concept is practical, providing the missions 
themselves are coordinated to the point of being clearly expressed, lacking in duplication and 
subject to a minimum of interpretation. When missions of this character are assigned, they 
breed the coordination necessary for effective support of forces. 

Examples 

This principle is well exemplified in the weapon system concept of management being 
practiced today within the military services. This concept of management on a total systems 
basis involves the superimposing of a systems management pattern over a basic functional 
management structure. The fundamental intent of this concept is to insure effective coordina- 
tion of the activities within each of the functional areas in order that the complexities and 
peculiar demands of modern weapons may be dealt with in a most effective manner. On this 
basis specific functions continue to be recognized as entities performing needed roles but their 
individual actions are continuously made known to each other. It is perhaps a first step toward 
integration (a separate principle) in which the functional areas may be administratively and 
physically organized as a single unit. 

Perhaps of even greater significance than coordination internally in logistics is the 
need for coordination between logistics and tactics. Innumerable examples spring to mind 
here—McArthur's campaign mentioned earlier, Gen. Patton's dashes through France and 
Germany, and of course, back in history, Napoleon's abortive campaign to Moscow and back 
provides a stark picture of the consequences of failure to achieve this coordination. 


PRINCIPLE OF SUPPORT INTEGRATION 


The effectiveness of logistics is in direct proportion to the degree 
of integration achieved among the various support functions. Forces 
require support as a unified effort in which each of their resource needs 
is provided through an integrated logistics system. Support which is 
disjointed obstructs force capabilities. 


Effective support of forces is not the product of piecemeal effort. It is the sum of 
mechanical means and human endeavor. It is born of an integrated system. 

Logistics support is both a part of a system, of military force, and a system itself. 
Its primary endeavors involve determining resource requirements to meet needs set by the 
military force; acquiring them; and distributing, maintaining, and ultimately disposing of them. 
Its functions involve research, development, supply, transportation, maintenance, procurement, 
and production of materiel; recruitment and training of personnel; design, construction, and 
maintenance of facilities; and planning of all resources preparation. Its effectiveness, factored 
by the economy with which that effectiveness was achieved, determines its relative efficiency. 

Efficiency, effectiveness, and economy are not forever synonymous, They coexist where 
related activities are meshed, where duplication is avoided, and where the process flow from 
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one activity to another and through the entire system approaches a simple pattern. Integration 
is the means; support effectiveness the end. 

The pattern of integration bears watching. That integration which is accomplished on 
the basis of better support of forces is desirable. That which is accomplished merely as an 
internal logistics matter, to produce greater economy, for example, must be judged first as to 
its effect on force support, not on economy gained. 

Examples 

A current example of the application of this principle exists in the mission of the Air Force 
Logistics Command of the USAF. Here, at least insofar as wholesale logistics is concerned, is 
an effort to consolidate—integrate if you will—the various functions involved in supporting an 
entire military service. While the functional areas themselves—procurement, supply, main- 
tenance, transportation, etc.—remain identifiable, they are grouped together under one com- 
mand structure in the interest of assuring maximum compatibility of their efforts. 

Still another example of this principle is evident in the still developing concept of 
"single managership" of selected resource areas. While the integration possible here may not 
reach the point of actual consolidation of entire organizations, it does nonetheless go well 
beyond the scope of the principle of coordination in its intent to preclude duplication and assure 
maximum use of available resources, 


PRINCIPLE OF EVOLUTION 


Logistics must be oriented to accept change through the process of 
evolution. Evolution permits alteration of logistics systems and tech- 
niques without chaos. Evolution in logistics is mandatory to accom- 
modate the changing nature of forces to be supported. 


In fulfillment of the objective, logistics must follow a definite pattern of development. 
This is so, because, while the end objective—support of operational forces—does not change in 
itself, the forces making up that objective are by their very nature dependent on change. Change, 
as a means of retaining combat effectiveness, is mandatory. But, the military force which 
creates the logistics objective does not change all at once. As portions of it are modified, 
modernized, or otherwise altered in character, other portions remain stable. The logistics 
objective must recognize the past, the present, and the future. It must accept change while 
providing for nonchange. 

The need for, and indeed the fundamental objective of, support of combat has not changed 
since man first waged battle. Only the techniques and organization for support have changed, 
and then only through evolution. How else could it be otherwise? To the state, war may be 
thought of in terms of limited, general, and total; to the individual (his life in jeopardy), war is 
always total. So, through the ages, military men have been unwilling to rashly discard proven 
means of supporting their needs in favor of new and perhaps somewhat radically different 
means. They have chosen, as a matter of self-preservation, to leave well enough alone for the 
moment and instead to add to and build on their existing support structure. 

There are obvious dangers in the process of such evolution. Primary among these is 
the hazard of complexity. Complexity is of itself not necessarily an evil; yet considering that 
objectives are met only through the activity of human endeavor, it becomes of paramount 
importance that the human endeavor not become submerged by the scope and density of 
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techniques through which the very endeavor is to be applied. The antithesis, simplicity, must 
be recognized as an objective—not a principle—within the principle of evolution. Simplicity 
cannot be taken as a principle of logistics because it is a relative thing—relative in the sense 
that it is the job to be done, not the means of doing it, which decides the degree of simplicity 
attainable. The more complex the job, the more stringent its demands, the less simple can be 
the support structure. The objective of simplicity involves a search merely for the simplest 
structure or system which will do the job effectively; the simplicity achieved through this 
search is a matter of degree. 

The principle of evolution, then, is to be pursued between the constraints of complexity 
and simplicity, with a definite leaning toward the latter. This leaning is present when we 
clearly recognize that evolution must be practiced not just as an adding to but also as a taking 
away—of outmoded techniques, of disproven practices. 

Examples 

The weapon system concept of management previously discussed may be cited here as 
an evolutionary technique designed to keep pace with the requirements of modern technology 
while still allowing for operation of methods designed for support of earlier technology. 

Support of naval forces while at sea offers another example of this principle. From the 
early days when the Fleet put into a limited number of major ports for replenishment, through 
the World War II era when advanced support bases throughout the Pacific were set up to reduce 
the "downtime" of ships, to the latest techniques of underway replenishment (including carrier 
onboard delivery), the whole series of refinements offer a striking example of the principle of 
evolution within the logistics process. 


PRINCIPLE OF DISPERSION 


Logistics support of forces at the time needed and in the quantity 
needed is best assured through dispersion of resources. Dispersion has 
the character of placing an appropriate share of resources directly in 
the hands of the force and of protecting the reserves from destruction 
through diffusion. 


It is not the total resources possessed by a nation which decide the combat; only those 
resources which are actually applied by military force are decisive once combat is actually 
undertaken. Dispersion places maximum resources, and potential, in the possession of the 
force before the force actually needs to apply them. This is the first and most significant rule 
of dispersion. 

There is obviously another reason for dispersion. This is to secure resources from 
destruction. Dispersion, in this sense, avoids concentration and vulnerability. But this aspect 
of dispersion can only be practiced to the extent that it does not detract from responsiveness 
to the needs of the force. 

This nature of dispersion actually supports the deterrent strategy and tactics of a 
nation, for through the protection thus offered to resources-in-being (the men, weapons, equip- 
ment, supplies, and facilities available to, but not necessarily in possession of, the military 
forces), the potential of a nation to wage combat is maintained. In this sense the resources are 
being applied in a war of men's minds. Again the interdependency of strategy, tactics, and 
logistics is demonstrated. In the case of actual combat the strategy and tactics employed 
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depend on dispersion of resources for success. The combat involving men's minds employs 
the potential represented in dispersed resources as a very real weapon; the resources are, 
in fact, the foundation for the strategy and tactics. 

Dispersion of resources, to support military forces, has been recognized since the days 
when troops "lived off the land."" A distribution, dispersion in effect, of resources was needed 
which would permit foraging parties to supply the maneuvering force without having to search 
great distances, thereby delaying the maneuver. This requirement to have resources dis- 
persed close to the operating units continues as a fundamental concept in modern military 
thinking. Similarly, dispersion of resources, as a means of reducing their vulnerability, has 
been practiced in every war. 

In a more recent context, the quantity and range of supplies needed to support military 
forces had expanded, during the past years, to a point almost negating mobility sought in 
military forces themselves. The alternative was to locate the logistics endeavor in separate, 
large depots, in positions most advantageous to support of forces—a form of dispersion. 

With the advent of modern weapons and weapon carriers, the vulnerability of these large 
depots, in their immobile state, gave further emphasis to dispersion as a basic principle of 
logistics. The current trend, again a modern application of this principle, is to place a major 
portion of logistics support in the hands of the user, thereby enabling him to fight with what he 
has in terms of dispersed logistics. 

Principles of war, specifically those of mobility and reserve, in fact dictate this 
principle of logistics. Mobility of forces is enhanced by dispersion of resources toward the 
user. Similarly, just as certain forces must be held in reserve as a backup for unforeseen 
possibilities during combat, so there must be an accompanying reserve of resources held ina 
dispersed situation. 

Application of this principle of dispersion must clearly recognize assurance of readiness 
and mobility of military force as the first requirement. Dispersion of resources to decrease 
their vulnerability can then only be accomplished to the extent that it does not detract from this 
condition. 

Examples 

Examples of this principle in actual use can be discerned in practically every military 
endeavor down through the ages. Admittedly, however, it has risen in significance in proportion 
to the increase in destructive potential available to the enemy. Evidences of the application of 
this principle up through World War I are still to be seen in the geographical distribution of 
military depots which can be traced, in part at least, to a desire to decrease the vulnerability 
of stores. More recently the military services have developed concepts (a) of "self-sufficiency" 
for their different forces and (b) of having the forces continually prepared to initiate and sustain 
combat on the basis of enough resources in their immediate possession; such a posture was 
made possible through dispersion of resources to the using units. 


PRINCIPLES OF FLEXIBILITY 


There must be a flexibility of logistics support to match the mobil- 
ity of forces. Responsiveness to the needs of operational forces is 
most readily assured through adaptability of logistics which is the basic 
measure of its flexibility. Since logistics exists solely to support 
forces, assurance of flexibility is of priority importance. 
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The mobility of forces and accompanying logistics must be enhanced by flexibility of 
the basic logistics system. The system must be readily adaptable to changes in strategy, 
tactics, and operational techniques of the supported forces. This flexibility must exist in all 
logistical endeavor and rightly constitutes a principle of logistics. 

Flexibility, or the capability to react rapidly and reliably to changing situations, is a 
mark of effective logistics. It offers to logistics characteristics similar to the mobility of 
combat forces. Flexibility is a mark of the energy of logistics. 

The support of forces—Logistics—must embody flexibility, since employment of the 
forces themselves involves constant flexibility. Flexibility of force—actually a product of the 
principles of economy of force, mobility, and coordination—is, under modern conditions, more 
and more determined by the flexibility of support. It is not enough merely to have the right 
resources; rather, the right resources must be at the right place at the right time. This isa 
foremost objective of logistics flexibility. It is an objective to be attained through responsive- 
ness, a condition of flexibility. 

For the opposite condition to exist, that is, for logistics to be inflexible, invites 
disaster. Inflexibility, in this respect, prohibits application of the war principles of initiative, 
concentration, mobility and coordination. A logistics system which is inflexible assumes, 
quite improperly, the decisive role in force employment, for the forces, instead of being able 
to operate wholly toward the objective, must conduct their missions within the restraints 
placed on them by logistics. While it is quite evident that modern logistics is assuming a 
significance far beyond that previously offered, we must not lose sight of the fact that logistics 
is never an end in itself. Logistics exists only to support forces. Victory in combat is 
achieved by the application of resources through combat forces, not by the preparation of 
resources through logistics. 

The final decision must always rest with those who apply the resources: the combat 
force. Strategic and tactical decisions—while most effectively made in the context of logistical 
capabilities— must, to the greatest extent possible, be decisions made with the operational 
objective clearly in mind. Adherence to the principle of logistical flexibility will enhance such 


decisions. Adherence is best assured by application of another principle: that of logistics 
control. 


Examples 

In his book, "Logistics in the National Defense,"' Rear Admiral Henry Eccles cites the 
situation in 1944 wherein General Eisenhower changed the tactical plans of the invasion of the 
Continent to permit General Bradley to turn east, away from his original objectives, the 
Brittany ports, in order to join with the British forces approaching Falaise. He quotes General 
Walter Bedell Smith: "But a hazard greater than the now thoroughly routed Germans was 
troubling us—supply. It is no great matter to change tactical plans in a hurry and send troops 
off in new directions. But adjusting plans to the altered tactical scheme is far more difficult. 
It involves relocating vast depots and stores of ammunition which must flow to the fighting 
troops in an uninterrupted stream. Our bombing of French rail centers which had contributed 
so heavily to victory in Normandy now returned to plague us. The railroads were practically 
unusable, We laid out two-lane, one-way motor routes across France over which the trucks 
roared day and night to keep the advance supplied. Even this was not fast enough for the racing 
armored spearheads. They got their supply almost entirely. by air." 
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Here was demonstrated a need for logistics flexibility in its truest sense. Here was 
demonstrated, in the ultimate success of this portion of the campaign, the validity of this 
principle of logistics. 


PRINCIPLE OF CONTROL 


The preparation of resources must be controlled to insure that the 
quality and quantity of resources provided the military forces and the 
time in which they are provided will permit successful fulfillment of 
the mission. Logistics control must always serve the purpose of the 
supported forces; that control which is applied for other purposes has 
lost the objective. 


Two aspects of control are present in logistics: (a) control of the specific techniques of 
logistics so as to assure effectiveness of the processes internal to the logistics system, which 
represents the "science" of logistics, and (b) control of the blending of all these processes into 
effective support of forces, which represents the "art" of logistics. Of these two control aspects, 
the first may be defined as a management control whose objective is perfection and integration of 
techniques; the second as a command control whose objective is transmission of the combined 
techniques into support of forces. These two controls are distinct, yet obviously complementary. 
Management control recognizes as its concern the effectiveness and economy of logistics com- 
ponents; command control as its concern the compatibility of logistics operations with support 
needs. Both obviously must operate toward the same objective: effective support of forces. 

The question of who should "commana" the logistic support effort is frequently debated 
to a length beyond reason. Experience has shown that the unit or area commander need not 
necessarily have command control of his logistic support to be assured of its effectiveness. In 
many respects, having such control proves more of a burden than a boon. 

Command control of the logistics effort can very soundly be lodged in a separate com- 
mander, provided the commander diligently observes all the other principles of logistics 
enunciated here. The fact is, command need not be all one way or the other. Command must 
be lodged where it can assure the most effective results. Its placement is a function of the 
particular conditions surrounding the situation in question. Whoever exercises the command 
control, and does so on the basis of the recognized principles, is reasonably sure of success. 

The question of command control is not so much who should control but, rather, how to 
exercise the control so as to create the proper atmosphere between supported forces and 
logistics activities. The support element must be controlled, in a command sense, in such a 
fashion as to assure adherence to the Basic Principles of Logistics. Specific emphasis must 
be directed to providing the operational force with the quality and quantity of resources needed, 
at the time needed, to successfully fulfill the operational mission. In this climate, the force 
commander must strive to know what he can expect of the logistics element, the logistics com- 
mander to know what is expected of him. To the extent the force tommander and logistics 
commander learn to work together, command control of logistics becomes most effective as 
virtually a joint effort. 

There is a point here to be noted. Logistics is not an endeavor assigned to and con- 
trolled by any one organization. It has no definite starting and stopping points in the military 
System. It exists and is performed at every echelon. That is to say, that up until the final 
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moment of application, resources are being prepared—an act of logistics. This preparation is 
accomplished at all levels, from the industrial structure of the nation, on down to the individual 
military combatant. Actually, there is a control aspect involved at every echelon. The prin- 
ciple of logistics control is the same at every echelon; the preparation of resources must be 
controlled in a manner designed to ensure that the quality and quantity of resources provided to 
the military force and the time in which they are provided are sufficient to permit successful 
fulfillment of the mission. 

The first aspect of control involves a management effort: a measurement of logistics 
to assess its efficiency in terms of economy and effectiveness. Measurement must be made 
internally and externally. Although internal measurement of the logistical activity may give 
considerable emphasis to economy, external measurement must emphasize effectiveness. Of 
these, the latter is by far the more important. External measurement, in terms of readiness of 
forces, is the only reason for logistics. 

Internal measurement is only a subsidiary designed to identify specific weaknesses 
within the logistics process which are detracting from the support of forces. Internal measure- 
ment is significant only as it relates to logistics effectiveness as indicated by external meas- 
urement. Control is properly applied when it measures, for the force commander and the 
logistics commander, logistics support in terms of readiness of force; it is inadequately applied 
when it indicates only the relative efficiency of internal logistics processes. 

Examples 

A recent illustration of this particular principle is offered by the decision of the Air 
Force in 1956 to reassign control of its oversea logistics support depots from the theater 
commanders to the Air Materiel Command. Obviously this action was to be taken with certain 


misgivings on the part of some; yet, the reassignment was made, and the support of forces of 
the theater commander did not suffer. 


From this application, the principle has subsequently been extended, through the "'self- 
sufficiency" concept, mentioned earlier, to place control of those resources which he needs 
immediately in the hands of the organizational commander. The intent of all these actions, in 
keeping with the fundamentals of the principle, is to achieve a maximum effectiveness. 

With respect to the internal aspects of logistics control, increasing emphasis within the 
logistics endeavor is being placed on such measurements as weapon readiness, equipment 
downtime, manpower skill levels, and the like, all of which are representative of the indications 
of support effectiveness. Perhaps lesser emphasis is being placed on such measures as 
(a) internal depot processing time, (b) quantities of items requisitioned which were not supplied 
within time limits, etc.; these are measurements of the internal efficiency of logistics. 


PRINCIPLE OF ECONOMY OF RESOURCES 


The essence of successful logistics is to do more with less— 
through economy of resources. Economy of resources seeks to avoid 
depletion while assuring that those resources authorized supported 
forces are readily available. Economy of resources is achieved not 
only through an exchange of quality for quantity but is predicated on a 
continual striving for most effective management of resources. 
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There are eight primary resources with which military forces must secure their mis- 
sions. Of these, people, supplies, and facilities may be designated as basic. Services, trans- 
portation, and communications constitute functional resources, Money, time, and technology 
may be construed as determinant resources. 

Concerning basic and functional resources, the virtue of accomplishing a logistics task 
with the least quantity seems obvious. Yet there is danger in interchanging the words "economy" 
and "least."" 'Least" may be too little while ''economy" suggests providing no more than the 
minimum amount and degree of support needed by the military force to do its job effectively. 
"Least" implies, primarily, quantity; ''economy," on the other hand, involves a combination of 
quality and quantity. The nature of "economical logistics" is not necessarily that of least 
quantity but involves an input of quality so that mission accomplishment is in fact enhanced 
rather than jeopardized. 

In this context the principle of economy of resources stems from the principle of con- 
trol and feeds directly into the principle of dispersion and into the principle of mobility of 
forces. In both cases, reduced quantity, permitted by improved quality, contributes to the 
dispersion and mobility criteria. 

The fact that little or no observance of this principle is apparent in World War I or II 
does not negate it. This apparent lack of observance brings to attention the existence of the 
resources of money, time, and technology. There will be some who will disagree that money is 
a resource, more who will argue that time and technology cannot be viewed as resources. 
Obviously, no one of them can be physically applied by military forces in the sense that the 
other resources can. Still, the money available and used by the United States in World War I 
and Il, coupled with the time available and used to prepare for combat and the technology of 
materiel at the time, tends to afford these the sense of resources. On such a basis, money, 
time, and technology are to be thought of as "determinant resources" in that they determine to 
a large extent the quantity and quality and proportions of the basic and functional resources 
which will be available to military forces. The amount of time available, coupled with the 
amount of money, will determine whether more or fewer people, more or fewer supplies, more 
or fewer facilities, or any combination, thereof, can and should be provided. The technology 
existing at the moment will determine the quality in which the basic and functional resources 
are made available. 

Whereas economy of the resources of people, supplies, facilities, and money has long 
been advocated, there seems to have been too little appreciation for the need for economy of 
time. Here again, it is the advent of modern technology which has prompted greater recogni- 
tion of the time element: time to keep militarily ahead of, or even up with, capabilities of 
potential enemies, time to prepare the civilian population for what may happen in the event of 
war, and time to reconstitute resources in case war comes. The question is not how much 
time have we, but rather how well—a quality aspect—do we use the time available tous? We 
must learn to manage time to assure an economy of this resource. For quantity of time, here- 
tofore available to us, we must substitute quality in our use of time. 

Examples 


Application, though perhaps not on an entirely willing basis, of this principle can be 
found in the breakout and pursuit phase of the Continental Invasion in 1944. In August, for 
example, the materiel staff at Supreme Allied Headquarters deemed it inconceivable that suf- 
ficient support could be provided to maintain the advance on the scale and at the speed with 
which it had thus far been carried. Through the most economical allocation of materiel 
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resources, it was estimated that only 11 U. S. divisions could continue the drive, and even then 
levels of maintenance would have to be drastically reduced—to 162 tons per day for infantry 
divisions and 247 tons for armored divisions fully engaged, for example. Advancing divisions 
were to take with them none of their heavy artillery and only 50 percent of their medium 
artillery, and ammunition expenditures were to be limited to one third of normal. That combat 
operations were able to continue, and even be carried well beyond those originally planned, 
offers a marked justification for the validity of this principle. 

The aspect of improved quality as a means of lessening quantity has been aptly apparent 
in the technological advances in weapons and support techniques over the years. As the 
weapons—the atomic capability and the missiles offer striking examples—and specific tech- 
niques of support—automatic requisitioning procedures, electronic data processing equipment, 
and the like—have improved the quality of our effort, it has been proved possible to reduce the 
quantity element of the military inventory. 


PRINCIPLE OF PERSONNEL PARTICIPATION 


Logistics is forever the product of human endeavor. The flexibil- 
ity, the integration, control, and coordination of logistics are more the 
result of human initiative, ingenuity, and integrity than of any other 
factor. To the extent that people participate enthusiastically and 
effectively, the logistics process is effective. 


Logistics, for all its mechanical and materialistic aspects, is fundamentally the product 
of people. People are both the cause of and the solution to problems in the logistics field. In 
terms of the quality/quantity aspects of economy of resources, they present somewhat of a 
dilemma. There is no substitution for quality of people, and yet too frequently the mistake is 
made of trying to make up for less quality through more quantity. This severely misshapes the 
situation and causes people to keep busy looking for things to keep themselves busy. The result 
is more complexicy, less effectiveness, 

The proper approach is for more quality and correspondingly less quantity. There is 
needed in the field of human endeavor, as in the field of mechanical equipment, an upgrading 
of technology—thought technology. People must be given the opportunity, and encouraged, to 
increase their powers of observation, of understanding, and of objective reasoning. In logistics, 
as elsewhere, the human mind is readily capable of translating ideas into machinery; it is not 
yet so capable of translating the potential of this machinery into practical work. When technology 
of human thought has reached the level of machine technology (and it should surpass it), the 
participation of personnel as the motivating force of logistics will be truly effective. There is 
an urgent need for professionalism in the field of logistics. 

Aside from initiative and ingenuity, the quality of individual integrity is the insurance 
factor in logistics. Too often integrity is confused with strict, and rather blind, adherence to 
the letter of the military regulation. The integrity needed in logistics is that of the individual 
who conscientiously applies himself to the most effective performance of his particular function, 
with due regard to moral and legal restraints placed upon it. He does not do things a certain 
way merely because they have always been done that way or because the regulation says they 
will be done that way. He does things a certain way because that way is best for effective sup- 
port of the forces. Certainly, if this way runs contrary to established procedures, he questions 
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the procedure, for this is an element of integrity. Integrity of the individual is the moral fiber 


of logistics. 

Despite human shortcomings, people remain—in many situations—the saving grace of 
logistics. They provide the "extra" which is needed to overcome deficiencies or to make up 
for mechanical failures. Through ingenuity and improvisations, the system is made to work, 
when from all outward appearances, it is incapable of working. 

The logistics commander who recognizes the potential of his people, who effectively 
guides this potential and seeks to increase it even more, is the commander most assured of 


success. 

Examples 

Situations wherein this principle is evident range the spectrum of the logistics process. 
From the relative simplicity of the employee suggestion program through the complexity and 
gravity of supporting a combat force under the usual, near chaotic conditions of the initial 
stages of attack, the participation of the individual is the dominant factor. The reader has only 
to recall his own innumerable experiences, he has only to reflect on the support requirements 
as they are known today and forecast for tomorrow, to appreciate the foundation for this 
principle. The logistics process is a creation of man's mind, its operation a result of man's 
ability, its success or failure a product of man's ingenuity. 











CONSTRAINT QUALIFICATIONS IN MAXIMIZATION PROBLEMS*! 


Kenneth J. Arrow, Leonid Hurwicz, and Hirofumi Uzawa 





Many problems arising in logistics and in the application of mathematics 
to industrial planning are in the form of constrained maximizations with 
nonlinear maximands or constraint functions or both. Thus a depot 
facing random demands for several items may wish to place orders for 
each in such a way as to maximize the expected number of demantds 
which are fulfilled; the total of orders placed is limited by a budget 
constraint. In this case, the maximand is certainly nonlinear. The 
constraint would also be nonlinear if, for example, the marginal cost of 
storage of the goods were increasing. Practical methods for solving 
such problems in nonlinear programming almost invariably depends on 
some use of Lagrange multipliers, either by direct solution of the 
resulting system of equations or by a gradient method of successive 
approximations (see [5], Part II), This article discusses a part of the 
sufficient conditions for the validity of the multiplier method. 














INTRODUCTION 

This article covers an examination of the interrelationships of the additional assump- 
tions under which the following two propositions, both extensions of the classical Lagrange 
multiplier method ([1], p. 153), are valid. 


Quasi-Saddle- Point Condition 
If X maximizes f(x), subject to the constraints g(x) 2 0, and f(x) and g(x) are dif- 
ferentiable, then there exists y 2.0, such that f, + y 8, = 0, y g(x) = 0.! 





*Manuscript received January 28, 1959. 

This work was supportedin part bythe Office of Naval Research (Task NR-047-004) at Stanford 
University. Hurwicz's participation was made possible by a Rockefeller Foundation grant to 
Stanford for mathematical research in the social sciences. 

Ix is a column vector with components x],...»5 Xp, y a row vector with m components, f(x) 
is a real-valued function of x, g(x) a column vector with real-valued components gJ (x), 


G= 1... 9k f, a row vector with components fx; = 0£/0x;, 8x a matrix with components 


gJ = 0 gj/a x; where i varies over columns and j over rows. Bars over f and g or their 
a 


derivatives denote evaluation at x. If v is a vector, v 20 meansthat each component of v is 
nonnegative; v > 0 means that each component is positive. 

A function f(x) of a vector variable is said to be differentiable at x, if there is a row 
vector a such that 


lim [f(x + h) - f(x) - ah]/|h] = 0. 
h—0O 
(Continued) 
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Saddle-Point Criterion 

If f(x) and g(x) are concave, then a necessary and sufficient condition that x maximize 
f(x), subject to the constraints, g(x) 20 is that there exist a y 20, such that (x, y) isa 
saddle-point subject to y 20, of the Lagrangian function f(x) + y g(x). 

It is known that neither proposition is valid without additional assumptions.* Kuhn and 


Tucker [11] showed that both propositions are valid if f and g are differentiable and the fol- 
lowing condition is satisfied. 





Constraint Qualification K T° 

For all X in the constraint set C (defined by the conditions g(x) 2 0) and all X such 
that ex - (X - X) 2.0 for each component k of g(x) for which eX (x) = 0, there exists a dif- 
ferentiable vector-valued function (0) such that ¥(0) = x, ¥(@) belongs to C for all positive 
@ sufficiently small, and W'(0) = x - x. 

Also in this article, further results on the subject are discussed and simplified proofs 
are given. First, Constraint Qualification KT are slightly weakened so that the meaning of 
the qualification becomes more transparent. Theorem 1 shows that the Lagrangian method can 
be applied to those constrained maxima for which the weaker version of the Constraint Quali- 
fication is satisfied. Next this article shows that the Constraint Qualification in the present 
formulation is the weakest requirement for the Lagrange method to be applicable; namely, in 
Theorem 2 below, it is proved that if the Lagrange method is justified for all differentiable 
maximands (or even all linear maximands), then the constraint function satisfies the Constraint 
Qualification provided the constraint set is convex. 

The direct verification of the Constraint Qualification in specific cases is difficult, 
and it is useful to find simpler hypotheses which imply it. Several apparently new conditions 
implying the Constraint Qualification and, therefore, the validity of the Saddle-Point Criterion 
and the Quasi-Saddle-Point Condition are proved in a later section (A Sufficient Condition for 
the Constraint Qualification). Note that for differentiable functions f(x) and g(x), the Quasi- 
Saddle-Point Condition implies the Saddle-Point Criterion. For if f(x) and g(x) are concave, 
the Lagrangian function, f(x) + y g(x), is concave in x for any given y 2 0; if X maximizes 





1( Continued) io 
If this condition holds, then the partial derivatives of f(x) all exist at x and f, = a, but the 
condition of differentiability at a point is stronger than the existence of partial derivatives. 

The function f is said to be concave, if f[tx + (1 - t) x] 2 tf(x) + (1 - t) f(x) for any pair 
of vectors, x, x, and any real number t, 0 Stl. 

A vector function f is said to be differentiable (or concave), if each component is. 

A saddle-point, subject to y 20 of a function L(x, y) of two vectors, is defined by the 
properties that x maximizes L(x, y) and y minimizes L(x, y) for nonnegative y. The con- 
cept is due to Kuhn and Tucker [11]; see their definition of the Saddle Value Problem on p. 482 
of Ref. [11]. A quasi-saddle-point is a point where the first-order (necessary) conditions for 
a saddle-point are satisfied. 

In this article, the variables are not necessarily restricted to be nonnegative. Any such 
restrictions are therefore assumed tobe included amongthe conditions g(x) 20. The formula- 
tion of the following conditions therefore differs in detail but rfot in essence from that of Kuhn 
and Tucker [11]. 

2See, for example, Courant, [7], pp. 189-190 and 192-93. For the case of inequalities, the fol- 
lowing example is dueto Slater [12]: f(x) = x, g(x) = -(1 - x)*;here both f and g are concave 
and differentiable, yet there is no saddle-point and hence no quasi-saddle-point. See also the 
example of Kuhn and Tucker, Ref. [11], pp. 483-84. 

3See Ref. [11], p. 483. For a’corresponding condition in the context of equality constraints, see 
Bliss, Ref. [6], p. 210, conclusion of Lemma 76.1. 
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f(x) subject to g(x) 2 0, the Quasi-Saddle-Point Condition implies that the Lagrangian function 
has a zero derivative at x and, therefore, as a concave function, must have a maximum there. 
Since the Lagrangian is linear in y and g(x)2 0, it is obvious that y g(x) = 0 implies that y 
minimizes f(x) + y g(x) subject to y2 0. The converse part of the Saddle- Point Criterion, 
that if (x, y) is a saddle-point of the Lagrangian for some y, subject to y 2 0, then x isa 
constrained maximum of f(x) subject to g(x) 2 0, holds without any assumptions on g(x). 

(See Kuhn and Tucker, [11], Theorem 2.) 

Still another constraint qualification has been given by Hurwicz ([5], Chapter 4, Section 
V.3.3.2). In another section (Equivalence of Constraint Qualifications KT and H) it is shown 
to be equivalent to Constraint Qualification KT, at least for finite-dimensional spaces. 

To state these results more precisely and to relate them to other work, we will intro- 
duce some notation and definitions. In the first place, we define the constraint set, 


(1) C = {x: g(x) = oO}. 


In the second place, we will denote x - x by &. Finally, to simplify the statement of 
the conclusion of Constraint Qualification KT which will appear frequently in the following 
discussion, we will find it convenient to introduce the following definitions: 

Definition 1. A contained path (with origin x and direction é) is an n-vector-valued 
function (6) of a real variable which satisfies: 





(2) (6) is defined for all 0 = 6 £0 for some 9 >0; 
(3) W(0) =x, W(6)eC for all OS 059; 
(4) (0) has a right-hand derivative at @ = 0 such that w'(0) = ¢.. 


Definition 2. An n-vector § such that there is a contained path with origin x and 
direction ~ will be referred to as an attainable direction at x. The set of attainable directions 
(at any given X) will be denoted by A. 

The set of indices {1, ... , m} is divided into two parts, E and F. E is the set of all 
indices effective at x, namely, 





(5) E = 


and F is the set of all indices ineffective at x, namely, 


(6) F = {j: g(x) > oO}. 


Definition 3. An n-vector ~ is termed a locally constrained direction if = é 20 
(i.e., z) — 20 for all j€E). The set of locally constrained directions will be denoted by L 





With these definitions the Kuhn-Tucker Constraint Qualification can be written. 


Constraint Qualification K T 
Every locally constrained direction is attainable, i.e., LC A. 
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(Since the definition of a contained path does not require it to be differentiable through- 
out, this formulation is apparently weaker than Kuhn and Tucker's [11]. We do not know if the 
weakening is more than apparent when g(x) is differentiable.) 

We now observe that the set L of locally constrained directions is a closed convex cone, 
(By a cone is meant a set which, if it contains any point x, also contains x for every scalar 
420.) The set A of attainable directions is a cone but is not necessarily convex. 

Definition 4, Let W be the closure of the convex cone spanned by A, the set of attain- 
able directions (i.e., the smallest closed convex cone containing A). The elements of W will 
be termed weakly attainable directions. 

A weakly attainable direction is, then, the limit of a sequence of nonnegative linear 
combinations of attainable directions. We now introduce a weaker constraint qualification: 





Constraint Qualification W 

Every locally constrained direction is weakly attainable, i.e., LC w.t 

It will be shown (Theorem 1) that Constraint Qualification W is sufficient for the 
validity of the Quasi-Saddle-Point Condition and therefore for that of the Saddle-Point Criterion, 





4To see that A is not necessarily convex and Constraint Qualification W is truly weaker than 
Constraint Qualification K T, consider the constraints x; 20, x2 20, -X]X2 20. The constraint 
set C consists of the origin and the two positive half-axes. If ®¥ is taken to be the origin, the 
set of attainable directions A is the same as C. The set of weakly attainable directions, W, 
is the convex cone spanned by this set; that is, the nonnegative quadrant. All constraints are 
effective, and 


—e = 
ben = 0 S. bs 
0 0 


so that the set of locally constrained directions L is defined by gi = 0, gf, 20, and thus is 
again the nonnegative orthant. L is therefore contained in W but not in A. 

In general,if A were a closed set, thenthe convex cone spanned byit wouldalso be closed, 
and the words, ''the closure of,"' in Definition 4 could be deleted without loss of generality. It 
can be shown fairly easilythat itis closed for convex constraint sets. The referee has supplied 
the following example, which shows that A is not necessarily closed for differentiable g (x) 
in general: 

g (x, 0) - x* ‘ 


g (x, = -r@ {@4 sin® (1/6) + [max (r - sec @ tan 6 ,0)]*} for y>0, 


r = (x* + y2)!/2, 0 s @ = arc tan (y/x) S$ 7/2, 


and the domain of definition is the nonnegative quadrant. 
Since g(x, y) = 0 everywhere in the domain of definition, the constraint set 


Cs { (x, y): g(x, y) 2 o} is the union of the segments, 
{ (x, y): y = x arc tan (l/n7), x20, x* y}, 


for all positive integers n. Thus, the attainable directions from (0, 0) include (1, arc tan (1/n7)) 
for all n, but not (1, 0), so that A is not closed. 
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We now turn to some other conditions in the literature which have been found to be suf- 
ficient for the Quasi-Saddle-Point Condition. In the usual treatment of the Lagrange multiplier 
method in the case of equality constraints (Courant [7], p. 198), it is required that the matrix 

Z, of the (partial) derivatives of the constraint functions with respect to the variables have a 
rank equal to the number of constraints. This condition has been extended to the case of 
inequalities in [3], p. 8. We write, in the present notation, the 


Nondegeneracy Condition 

The rank of e* equals the number of effective constraints. 

(The condition given in [3] is actually slightly weaker.) It was shown in [3], Appendix I, 
that the Nondegeneracy Condition implies Constraint Qualification KT. Ina later section (A 
Sufficient Condition for the Constraint Qualification), we will deduce the Nondegeneracy Con- 
dition from a more general sufficient condition for Constraint Qualification W. 

In concave programming—that is, where the functions f(x) and g(x) are assumed to be 
concave—there are theorems which state conditions under which the Saddle-Point Criterion is 
valid but which do not involve Constraint Qualification KT. The simplest case is that of linear 
programming where the functions f(x) and g(x) are assumed to be linear. In this case, the 
Saddle- Point Criterion always holds, as is well known(({10], Chapters XIX and XX). Since the 
constraint qualifications are always assumptions about the constraint functions only and do not 
involve f(x), the question arises whether or not it is the linearity of g(x) that is vital. The 
answer is in the affirmative, in the sense that the linearity of g(x) is sufficient for the Quasi- 
Saddle-Point Condition; see Corollary 2 of Theorem 3. 

Another constraint qualification for concave programming has been proposed by Slater 
[12]: 


Constraint Qualification S 

The function g(x) is concave and, for some x’, g(x’) > 0. 

Slater showed that if the function f(x) is concave, his constraint qualification implied 
the Saddle-Point Criterion. A simplified proof was given by Uzawa ([5], Chapter 3). Slater's 
theorem was extended to more general spaces by Hurwicz ([5], Chapter 4, Theorem V.3.1). 
(Slater's assumption of continuity is dispensable.) 

Karlin ({9], Chapter 7, Theorem 7.1.1) suggested still another constraint qualification: 


Constraint Qualification K 

The function g(x) is concave and, for every y =0, there is an x such that y g(x) > 0. 

This condition is clearly implied by Constraint Qualification S; Hurwicz and Uzawa ([5], 
Chapter 5) have shown that in spaces of considerable generality the two conditions are in fact 
equivalent. 

It is natural to investigate the relation between Constraint Qualification S (or K) and 
the more general Constraint Qualifications, such as KT or W. Obviously, conditions S or K 
do not require differentiability of g(x), so that they cannot be completely subsumed under con- 
ditions KT or W, which do. We may, however, ask whether or not the former do imply the 
latter conditions under the additional assumption that g(x) is differentiable. In the section 
devoted to A Sufficient Condition for the Constraint Qualification, it is shown that such is 
indeed the case; in fact, a single condition (Theorem 3) is given under which all previous con- 
ditions can be subsumed as special cases, as well as some conditions not previously given in 
the literature. 
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A still weaker version of this condition is presented in Theorem 4. In one of its corol- 
laries, the hypothesis refers to functions which are simultaneously concave and quasi-convex, 
A characterization of such functions and also of functions which are simultaneously quasi- 
concave and quasi-convex is presented in the section, A Characterization of Functions Simul- 
taneously Concave and Quasi-Convex; this result may be of interest in other contexts. 

In [5], Chapter 4, Section V.3.3.2, Hurwicz introduces the following constraint quali- 
fication: 


Constraint Qualification H 

For all x€C and all & such that %. —€ +g(x)20, é is attainable. 

Since Constraint Qualification H does not identify separate coordinates, it is meaningful 
in all linear topological spaces in which a differentiation operation can be defined.’ It has been 
shown ([5], Chapter 4, Theorems V.3.3.2 and V.3.3.3) to be a sufficient condition for both the 
Saddle-Point Criterion and the Quasi-Saddle-Point Condition in spaces of considerable gener- 
ality. We will show in the section, Equivalence of Constraint Qualifications KT and H, that in 
finite-dimensional spaces, Constraint Qualifications KT and H are equivalent. 


PRELIMINARY LEMMAS AND REMARKS 
LEMMA 1: Every weakly attainable direction is locally constrained. 


PROOF: Let é be an attainable direction, (6) a contained path with origin x and 
direction — . Then for some @ > 0 and every jeE, 


gl [w(o)}=0, and gi fWe)})=0 (0<0< 4%); 


hence, for @= 0, dg! [W(6)//de = g) y' (0) = gi & 20 for every jeE, sothat é is locally 
constrained.© That is, A is included in L. Since L is a convex cone, the convex cone spanned 
by A must also be included in L; and since L is closed, W, the closure of the convex cone 
spanned by A, must also be included in L. (Q.E.D.) 


Let us define, for x€C, 
(7) K = closure of the set {A(x - xX): 4 20, xeC}. 


K is the union of all half-lines from x through elements of C, together with the boundary of 
the union. K is clearly a closed cone. That the set in braces is not necessarily closed is 


illustrated by the case g(x, y) = y - x’, x = (0, 0), where the x-axis belongs to K but not to 
the set in braces. 


LEMMA 2: If the constraint set C is convex, then K is a closed convex cone, and 
K CW. 


PROOF: The convexity of K follows immediately from that of C. 





It is possible, although we have not investigated this point, that Constraint Qualification KT 
can be extended in a natural form to infinite-dimensional (function) spaces. 

6Note that from the differentiability of gJ(x) and from that of Y(0@) at @ = 0, the chain rule for 
differentiation is valid. Theorem 6.14, p. 113 in [1] may be easily extended to the present case 
of one-sided differentiation. 





CONSTRAINED MAXIMIZATIONS 
If xe C, then by the convexity of the set C, 
x + O(x-x)eC forall Of @21. 
Hence, x - X is attainable and therefore weakly attainable. Since W is a cone, A(x - x) €W 
for all 420. (Q.E.D.) 
Let B be any set of vectors. The negative polar cone, to be denoted by B', is defined by 
B' = {u: ux £0 forall xeB}. 
We have (Fenchel [8], pp. 8-10), 
(8) B' is a closed convex cone; 
(9) B, © By implies that BY > By : 
(10) if B is a closed convex cone, B" = B. 
LAGRANGE REGULARITY AND THE CONSTRAINT QUALIFICATION 


Definition 5. An m-vector-valued function g(x) will be termed Lagrange regular if, 
for any differentiable function f(x), the Quasi-Saddle-Point Condition holds. 





LEMMA 3: If X maximizes f(x) subject to x€C, then 
f,eW' 


where W' is the negative polar cone of W. 


PROOF: Let W(@) be a contained path with origin x and direction é; then 


f[w()] s f[w(0)] = f(x) for all 05026. 


1, § f, ¥' (0) <0, 


for any — in A and, by continuity and convexity, for any ~ €W (Definitions 2 and 4). (Q.E.D.) 


THEOREM 1: If g(x) satisfies Constraint Qualification W, then g(x) is Lagrange 
regular. 

PROOF: Let f(x) be a differentiable function and xX maximize f(x) subject to x€C. 
Then, by Lemma 3, few". On the other hand, Constraint Qualification W states that L C W 
and therefore implies, from (9), that 


Hence, we have 
(11) 
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If B is the closed convex cone consisting of all vectors y= (-g ) with y= 


Definition 3 and (11) show that 4 € B", and by (10), 


20, then 


7 _sE-E <=E . 
4, , Se for some y 20. 


Define 


F 
y 


(y= FF) with pF <0. 





Then : + y 8, = 0, yg(x) = y= ge (x) = 0, from (5), so that the Quasi-Saddle- Point Con- 
dition is satisfied. (Q.E.D.) 

Theorem 1 is the basic necessity theorem for nonlinear programming ([11], Theorem 1) 
extended to the weaker Constraint Qualification W of this paper. 


THEOREM 2: If g(x) is Lagrange regular and if the constraint set C defined by it is 
a convex set, then g(x) satisfies the Constraint Qualification W. 


PROOF: It will be shown first that 
(12) ret . 
Let aeK'; then from (7), for A=1, 


(13) a(x-x) 20 forall xeC. 


Then X maximizes the function f(x) = ax subject to xeC. By the Lagrange regularity of g(x), 
there is an m-vector y such that 


(14) a+ye,=0; y20, 
and 
(15) ye(x)=0. 


The conditions (14) and (15) imply that y Fo and, thus, that 
(16) a+yegeao, yF zo. 


The condition (16) implies that 








aé =0 forall ~ such that gre 


aeL. 









Hence, we have the relation (12). Then, by Eqs. (9) and (10), K > L. Applying Lemma 2, 


WO>KD2L.  (Q.E.D.) 
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A SUFFICIENT CONDITION FOR THE 





” CONSTRAINT QUALIFICATION 
THEOREM 3: If E' is the set of effective constraints which are mee functions, E" 
is the set of all other effective constraints, and if there exists : such that ge :* 20, 
gee” ~ > 0, then Constraint Qualification W holds. 
PROOF: Let — be any element of L, @ any positive real number, and 
w(o)=x%+(— +a € JO, for 6 20. 
x 
Con- We will show that (6) is a contained path for @ sufficiently small and, hence, +a@ é is 
attainable. If we let a approach zero, it will follow from Definition 4 that — belongs to W, in 
rem 1) fact to the closure of A, so that we will have shown that L C W, which is Constraint Qualifi- 
cation W. 
For any j€E, at 6 = 0, 
it is 
P _s e - a a 
dg’ [v(0)]/de=Ei(E +a E)= gL Eragi se zagle, 
‘ ' 
since e) — = 0 by definition of L. If j€E , it follows from the hypothesis of the theorem that 
dg! [W(6)]/d@ 20 at 6 =0, 
which, for a convex function, implies that gl [Y(@)] has its minimum for @ 2 0 at 6 = 0. 
If jeE , then 
f g(x), dg! [W(e)]/a4e > 0 at 6=0, 


so that gl [¥(@)] has a local right-hand minimum at 6 = 0. It follows that 
E . E — ae 
g[W(6)] 2 ¢°[W(0)] = g(x) =0 for @ sufficiently small. 


Since, by definition, eF [w(0)] > 6, eF [w(6)] 2 0 for 6 sufficiently small, so that y(@)eC 
for 6 sufficiently small, from which the theorem follows, as has previously been shown. 


COROLLARY 1: If g(x) is convex, then g(x) is Lagrange regular. 


" a 
PROOF: In this case, Eis the null set, and it suffices to set € = 0. The conclusion 
follows from Theorem 1. 


As a special case, we state 





COROLLARY 2: If g(x) is linear, it ” Lagrange regular. 


COROLLARY 3: If g(x) is concave, E the set of effective constraints i. are 
co E' the set of all other effective constraints, and there exists x’ such that ge ‘(x )20, 
BY (x ) > 0, then g(x) is Lagrange regular. 
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PROOF: For any concave function, 


Ei (x" - x) 2 h(x") - h(x) = ite") for jek. 


Since the only functions which are both concave and convex are linear, the results follow from 
Theorem 3. 
The following special case is precisely Constraint Qualification S. 


* 
COROLLARY 4: If g(x) is concave and g(x ) > 0 for some x’, then g(x) is 
Lagrange regular. 


Corollary 4 was originally stated as Theorem 3 in Ref. [4]. The following corollary 
generalizes Theorem 2 in Ref. [2] which, in turn, generalized Corollary 4. 


COROLLARY 5: If the constraint set C is convex and possesses an interior, and 

et ~ 0 for each jeE, then g(x) is Lagrange regular. 
PROOF: By Lemma 2, x - x is weakly attainable for all x€C and therefore belongs to 

L by Lemma 1. Since C possesses an interior, L must possess one also and therefore has 
the full dimensionality of the entire space. If, for some jcE, 3 € =0 for all é in L, it 
follows that ze) € = 0 for all € , which means that ze = 0, contrary to hypothesis. Hence, for 
each j€E, there exists ! €L such that a) tJ) 20; since ee 2 0 for all € in L by defini- 
tion, we must have 


giel>o, gle*20 tor j,keE. 


If we let 


we see that 


and the conclusion follows trivially from Theorem 3. 


* 
REMARK: To establish that C has an interior, it is sufficient that g(x ) > 0 for some 
* 
x ; to establish that C is convex, it is sufficient that g(x) be quasi-concave.? 


We can also relate the Nondegeneracy Condition to this analysis. 


COROLLARY 6: If the rank of e® equals the number of effective constraints, then 
g(x) is Lagrange regular. 


PROOF: Let u be an arbitrary positive column vector; then from the hypothesis, we 
* ee * 
can find — — such that g® z = YS. 





A function f(x) is quasi-concave if, for all c, the set {x: f(x) 2 c} is convex. A function is 
quasi-convex if it is the negative of a quasi-concave function, so that the sets {x: f(x) = c} are 
all convex. 



















from 


ngs to 
has 

t 

e, for 
lefini- 


> some 


we 


ion is 


share 














CONSTRAINED MAXIMIZATIONS 185 





If we reconsider the proof of Theorem 3, we see that the convexity assumption for the 
elements of E' is only needed to insure that g! fy (6)] has a minimum at @ = 0 for jeE'. For 
the purposes of the proof, a local minimum is sufficient. It would therefore suffice to define 
E' as the set of effective constraints which are locally convex (i.e., which are convex over 
some neighborhood of x), which, for example, would be implied by well-known conditions on 
the matrix of second partial derivatives of the g! (x) evaluated at x = x. It is clear that the 
larger E' is, the weaker the hypothesis of the theorem. 

A somewhat different weakening of Theorem 3 is suggested by observing that if 6! is 


the ja unit row vector, then 


(17) -gi +56) gE =0 for any jeE. 


Since 6! eg (0) = 0, we see that the Lagrangian conditions for a constrained maximum of 
- ol (x + —&) subject to ge é 20 are satisfied at & = 0. Suppose for the moment that the 
Lagrangian conditions (i.e., those in the Quasi-Saddle-Point Condition) were sufficient to 
insure a constrained maximum. Then gi (x + —) would _ a minimum at é = 0 forée L. 
Since ( +a. &*) 6€L for all 6 2 0, this implies that g! [w(e)] 2 gl [Y(0)]= 0, and so the 
argument of Theorem 3 is still valid if E' is defined to contain all effective constraints for 
which the Lagrangian conditions are sufficient to insure a constrained maximum for - g! (x). 
A set of hypotheses, under which the Lagrangian conditions are sufficient for a con- 
strained maximum, is presented in [2], Theorem 1. We state the relevant results as a lemma 
(the statement below is slightly more general in that the domain of definition is not restricted 
to the nonnegative orthant; the previous proof extends easily). 


LEMMA 4: Let f(x) be a differentiable quasi-concave function and g(x) an m-vector 
valued differentiable quasi-concave function defined over some convex domain D. Let x and 
y 20 satisfy the Lagrangian conditions, f, +y Te =0, y g(x) = 0, and let one of the following 
conditions be satisfied: 

(a) f, x > % x? for some x! eC, x7eD; 


(b) . ~ 0 and f(x) is twice differentiable in a neighborhood of x ; 


(c) f(x) is concave. 
Then x maximizes f(x) subject to the constraints g(x) 2 0 
Lemma 4 can be applied to the preceding argument, with x replaced by é , f(x) by 
-g)(x+ &), and g(x) by ge &. Since the latter is linear and therefore necessarily quasi- 
concave, the hypotheses of the lemma are equivalent to requiring that one of (a), (b) or (c) 


hold for the function - g!(X+ &). The application of (b) and (c) is straightforward. Condition 
(a) becomes 


oi gi <i 23 
(18) g, § se : 


for some ¢} in the constraint _ '. E; 20 and some ¢2 for which x + Eg is in the domain 
of definition of Pa (x). Since ¢ _§ z 0 for all — eL, while OeL, (18) is equivalent to 







ei (x-x) >0, 





186 K, J. ARROW, L. HURWICZ, AND H. UZAWA 


for some x in the domain of definition of Pal (x). We can thus state the following generalization 
of Theorem 3: 


THEOREM 4: Let E' be the set of effective constraints g (x) which are quasi-convex 
functions and which satisfy one of the following three conditions: 


(a) ray x> glx for some x for which gl (x) is defined; 
(b) g) ~ 0 and Pal (x) is twice differentiable in a neighborhood of x = x; 


(c) g! (x) is convex. 


If E" is the set of all other effective constraints and if there exists "e such that = é * 29 
g= &*>0, then Constraint Qualification W holds. 

From this theorem can be deduced generalizations of some of the corollaries to 
Theorem 3. Corresponding to Corollary 1 of the latter, we have 


’ 


COROLLARY 1. If g(x) is quasi-convex and each effective component satisfies one of 
the conditions (a), (b), or (c) of Theorem 4, then g(x) is Lagrange regular. 


Corollary 3 of Theorem 3 can be generalized to 


COROLLARY 2: Suppose g(x) is concave. Let E’ be the set of effective constraints 
which are also quasi-convex and which satisfy one of the following conditions: 


(a) gl (x) > 0 for some x; 


(b) g ~ 0 and gl (x) is twice differentiable in a neighborhood of x = x ; 


(c) e (x) is linear. 


' 
If E" is the set of all other effective constraints and if there exists e* such that ge (x*) = 0, 
gE (x*) > 0, then g(x) is Lagrange regular. 


PROOF: If (a) holds, then, since Pal (x) is concave, 
B(x - %) 2 P(x) - (x)= Pw >0, 


so that (a) in Theorem 4 holds. Since (b) is the same in the two statements and (c) here implies 
(c) in the statement of Theorem 4, E' as defined here satisfies the conditions of Theorem 4. 
The proof is then the same as that of Corollary 3, Theorem 3.8 


A CHARACTERIZATION OF FUNCTIONS SIMULTANEOUSLY 
CONCAVE AND QUASI-CONVEX 
To better understand the domain of applicability of Corollary 2 of Theorem 4, it is 
useful to give a characterization in Lemma 6 of functions which are simultaneously concave 
and quasi-convex. First, we characterize functions which are simultaneously quasi-concave 





Sit should be stated that the case of nonlinear equality constraints is not well handled by these 


theorems and corollaries; these all depend on the construction of a linear contained path, but 
for nonlinear equality constraints none may exist. A form of Corollary 6 to Theorem 3 does 
remain valid in this case ([3], Appendix 1). 
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and quasi-convex. In what follows, an indifference set is a set {x: f(x) = c}; a maximal 
(minimal) indifference set is the set on which f(x) attains its maximum (minimum). A set S 
will be said to be bounded by two noncrossing hyperplanes in D if there exist linear functions, 
L, (x), Lo (x), not identically constant in D, such that 


(19) S = {xeD: L,(x) 20, Lo(x) =O}, 
and 
(20) L, (x) <e. La (x) > 0 forno xeD. 


A diagram will show that (19) and (20) are algebraic transcriptions of the geometric concept 
they define. 


LEMMA 5: A function is both quasi-concave and quasi-convex over a convex domain D 
if and only if every indifference set not minimal or maximal is bounded by two noncrossing 
hyperplanes in D. 


PROOF: Without loss of generality, assume that D is the domain of definition of f(x). 
We first note that quasi-convexity implies that {x: f(x) < ct is convex for all c. Suppose 
f(x°) <¢, £(x}) a x? a convex combination of x° and x~. If c' = max (f(x°), £(x})], we 
have x°, x- belonging to the set {x: f(x) < c'}, which is convex by definition of quasi- 
convexity, and hence f(x) S$ c' < c. 

We suppose, without loss of generality, that the linear space is the smallest containing 
D. Consider any value, c, of f(x), which is neither the maximum nor the minimum. If the 
set {x: f{(x) 2 c} did not have the full dimensionality of the space, there would exist a linear 
function L(x) not identically zero, such that L(x) = 0 whenever f(x) 2 c. Choose y so that 
L(y) = 0, and let z= y- x. Then, for any x, L(x +tz) isa linear function of t which is 
nonzero for t = 1 and therefore takes on the value 0 for at most one value of t; hence we can 
find € arbitrarily small for which L(x +e€z) = 0, L(x - ez) =~ 0. Suppose x is an interior 
point of D. Choose € sufficiently small so that x + €zeD. Since L(x+ ez) = 0, f(x+€z) <¢; 
by the convexity of {x: f(x) < c}, {(x) < c. Since D is convex, it follows by continuity that 
f(x) = c for all x€D, so that c would be the maximum value of f(x), contrary to assumption. 

The set {x: f(x) 2 c} is convex and has the full dimensionality of the space; the set 
{x: f(x) < cl has been shown to be convex. Hence, there is a separating hyperplane, i.e., a 
linear function, not constant over the entire space and hence not over D, L, (x) such that 
L, (x) = 0 for f(x) < ¢, L, (x) 2 0 if f(x) 2 c ([8], Theorem 28, p. 48), If L, (x) = 0 when- 
ever f(x) 2 c, the set {x: f(x) 2 c} would lie in a hyperplane and therefore not have the full 
dimensionality of the space, a contradiction. Hence, L, (x°) > 0 for some x°® for which 
f(x°) = c. Suppose L, (x) = 0 for some x for which f(x) < c. Then for y a convex com- 
bination of x° and x sufficiently close to x, we have L, (y) > 0, f(y) < c, contrary to the 
separation result. Thus, 


L, (x) <0 if f(x) <c, L, (x) = 0 ff £&) 2c. 


Since the sets {x: f(x) < c} and {x: f(x) = c} together exhaust D, we have 
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(21) {x: f(x) < ec} 


{x €D: L, (x) < 0}, 


(22) {x: f(x) 2 c} = {xeD: L, (x) 2 Oo}. 


Similarly, from the quasi-concavity of f(x), we find there exists a linear function, Lo (x), 
such that 


(23) {x: £(x) 


lA 


c} = {xeD: L(x) = 0}, 
(24) {x: f(x) > c} = {xeD: Ly (x) > 0}. 


From (22) and (23), (19) holds for the set S={x: f(x) = c}. Since the set (24) is included in 
the set (22), (20) holds. 

We now prove the converse theorem. First, consider any c which lies strictly between 
the maximum and minimum values of f(x). By assumption, the set S={x: f(x) = c} satisfies 
(19) and (20) for some linear functions, L 1), L 9 (x). Let S={xeD: L 1 &) < 0}, 

S = {x €D: Lo (x) > O}. We will first oe. that the set {x: f(x) = c} is nll Since 
c< _— tx, {(x° ) - c > O for some x®, , 


nm (19) we must have f(x) - c «0 for all x¢S, andall x eS. Further, f(x) -c 
cannot assume both signs in S, for, by the convexity of the set and the continuity of f(x), we 
would have f(x) - c = 0 in that set for some x, which has just been shown impossible. 
Similarly, f(x) - c must have a single sign in S. 

Suppose f (x°) -c>O0 for some x° in S. We will show that it is impossible that 
{(x") - c > 0 for some x! in S. For then, f(x) - c > 0 for all x¢S and all xeS, while 
f(x) = c for all xe S, and therefore f(x) 2c for all x €D, contrary to assumption. 

Hence, if £(x°) >e for some x® in S, f(x) > c for all x eS, and only for such x. 
The set {x: f(x) = c} is then precisely the set{x¢D: L 1 (x) = 0}, which is certainly convex. 

Alternatively, we might have f (x°) -c>0 for nia x°«S. The argument is com- 
pletely parallel. 

We have shown that {x: f(x) = c} is convex for any c, neither maximal nor minimal. 
The proof that {x: f(x) 2 c} is convex for such c is completely parallel. There remain the 


cases where c = max f(x) or min f(x). In the first case, the set {x: f(x) < c} is the entire 
x x 


set D and is certainly convex. The set {x: f(x) 2 c} is the intersection of all the sets 
{x: f(x) 2 c'} for c’ < ¢; it is an intersection of convex sets and therefore convex. The 


case where c = min f(x) is handled by the same argument. 
x 


LEMMA 6: A concave function is also quasi-convex over a convex domain D if and 


only if every indifference set not maximal or minimal is the intersection of D with a hyper- 
plane. 


PROOF: Let S, as before, be an indifference set, {x: f(x) = c}. We consider five 
cases: 


ce) 


(a) For some x’, x! ¢ s; L 1 (x°) > 0, Le (x! ) < 0. From (19), Ly (x° )<s 


L, (x ) 2 0. In this case, let x" = “(x2 + xl) /2, which belongs to S by a” 
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Ly (x’) > 0, La(x’) < 0. 


For any x €D, let x(@) = (1 - 6) x + 6 x”, g(@) = f[x(6)]. Clearly, x(@) €S for @ in the 
neighborhood of 1, so that g(6) = c for all 6 in the neighborhood of 1, and g'(1) = 0. Since 
g(6) is concave, it has its maximum at 6 = 1, and, in particular, f(x) = g(0) = g(1) = c, so 
that c = max f(x). 

x 


(b) L, (x) = 0 for all x €S, La (x) = 0 for all x é€D for which L, (x) 
case, x €S, if and only if xe D, L, (x) = 0. (Q.E.D.) 


0. In this 


(c) La (x) = 0 for all x€S, L, (x) 2 0 for all x¢€ D for which La (x) = 0. In this 
case, S = {xe D: Ly (x) = 0}. 


(d) L, (x) = 0 for all xe S, L (x°) >0, L (x°) = 0, for some x°e€D. First suppose 
1 2 1 ° 
L, (x) < 0 for some xeéD. Then we could find a convex combination of x~ and x for which 
L, is negative, Lo positive, contrary to (20). 


(25) L, (x) 20 forall xeD. 


Since L, (x) is not identically 0 in D, we can find x! so that L, (x) > 0. If Ly (x) <0 
for some x € S, we can find a convex combination of x, x! for which L, is positive, Ly nega- 
tive, contrary to the assumption that L, is zero for all xe€S. 


(26) Lo(x) = 0 forall xeS. 


From (25) and (26), the hypotheses of (c) are satisfied. 


(e) Ly (x) = 0 for all xe€S, L, (x°) < 0, La (x°) = 0 for some x°€D. This case is 
completely parallel to (d). 

Thus, except in case (a), the indifference sets are the intersections of hyperplanes 
with D, 

The converse follows trivially from Lemma 5; by assumption, for every c not maximal 
or minimal, there exists a linear function L(x) such that 


{x: f(x) = c} = {xe D: L(x) = 0}, 
so that (19) and (20) are satisfied with L, (x) = Ly (x) = L(x). 


EQUIVALENCE OF CONSTRAINT QUALIFICATIONS KT AND H 


THEOREM 5: In finite-dimensional spaces, Constraint Qualifications KT and H are 
equivalent, 


PROOF: Clearly, the hypothesis of Constraint Qualification KT can be inferred from 
that of condition H by considering only those components j for which g) (x) = 0. Since the 
conclusions of the two Constraint Qualifications are the same, the KT condition implies con- 
dition H. 










(1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 
[9] 


(28) 


(29) 


K. J. ARROW, L. HURWICZ, AND H. UZAWA 


To establish the converse, suppose Constraint Qualification H and the hypothesis of 






















Constraint Qualification KT hold. Then for any « > 0, 


if g)(X) = 0, e-(c €) 20. 
On the other hand, clearly 


if (x) > 0, g(x) + gl -(e —) 20 for ¢ sufficiently small. 


Since x €C, gl (x) > 0 for all j. From (28) and (29), g(x) + g(e ) 2 0. By Constraint 
Qualification H, there exists a path y _ (9) such that y, (0) = x, Vv, (8) eC for @ sufficiently 
small, Vi, (0) = « &. If we now define (0) = V.(6/ £), we have a contained path at x in 
direction é. 
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A DIRECT METHOD FOR THE COMPUTATION 
OF WAITING TIMEs*t 


Lionello Lombardi 





A recurrence relation for the direct computation of the waiting times 
of an array of c customers served by s parallel channels on a first 
come first served basis, where the time of arrival and the holding time 
of each customer are known, is worked out. The number of operations 
required by this procedure is proportional to c+ s, while the method 
commonly used, which is due to F. Pollaczek, requires a number of 
operations proportional to c2- In(c): this improvement is obtained by 
means of the adoption of sorting by insertion and by saving certain 
intermediate results. 











STATEMENT OF THE PROBLEM 

Let us consider a set of s equivalent parallel channels, and a sequence {C3} of cus- 
tomers who utilize the facility offered by these channels; the customers are served on a first 
come first served basis. Each channel may not serve more than one customer at a time, and 
each customer is served by the first channel available to him. No customer requires the 
services more than once—i.e., if any of them would do this, this model would consider him as 
anew customer each time he requires the service. 

The customer C; arrives at the entry-point of the channels at the time A and the 
service he requires is carried out in a time Tj: The sequence {A} will not be decreasing; no 
other assumptions about the sequences {A,} and {T,} are necessary in what follows. 

We want to compute the time Wi that the customer C; has to wait before he is admitted 
into one of the channels for service, for the first c values of the index i. 


POLLACZEK'S PROCEDURE AND ITS HANDICAPS 
Let us first point out that 


(2.1) Win = 0» (m = 1, 2,..., 8); 
i.e., the first s customers do not have to wait. 

For i > s the computation of Wi may be carried out using a method first introduced 
by F. Pollaczek [3]. However, Pollaczek's procedure was not really designed for being 
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directly used for the computation of Wi; it was rather intended to be a step of a more general 
scheme. On the other hand, we are emphasizing the computational aspects of the procedure, 


Let us briefly describe it. If E i denotes the time at which the customer C; leaves the 
channel where he was served, the equation 


(2.2) E, = A, + W, + T; 


holds. 


In order to compute W,, (i> s), the following three steps are performed: 
(a) The differences pii) = E, - A;, (¢@ = 1, 2,...,i-1) are computed. 
(b) The i-1 numbers pf) are sorted by decreasing order. 


(c) The sequence Bi), where 
- ( . (i) 
Bi = max [ 0; Di 
is computed. {pt is also not increasing. Then 


BY) = w,. 


In order to evaluate the efficiency of this method, we may first point out that the steps 
(a) and (c) require, for each i,(i = s +1, s+ 2,..., cc), a number of operations proportional 
to i. The number of operations which are necessary in order to carry out step (b) depends on 
the sorting method adopted; at the very best it will be asymptotically proportional to i - In(i) 
(see Bibliography [1], [2)]). 

In order to obtain the sequence {W,}, (i= 1, 2,..., c), the steps (a), (b), (c) have to be 
repeated (c - s) times for increasing values of i: hence the number  4( c) of operations 
which are necessary for calculating {Ww} is approximated by the formula, 


(2.3) $1 (c) = c? (Ky + In(c) - Ky) + f(s), 


where K, and Ky are positive constants and f(s) is a function of s. For large values of c, 
?4 (c) may be approximated by means of the formula, 


(2.4) ,(c) z= K+ c? + In(c), 


where K is a positive constant. 

A second burden of Pollaczek's procedure is the need of always having direct access 
to all the numbers Ai, Tj, Wi; this implies that 3(c - s) storage locations have to be 
reserved for this purpose; storage space is thus wasted, and if the size of the random access 
storage of the computer used is small compared to c, time-wasting dynamic transmission of 
information becomes necessary. 

A simple scheme that we might consider in order to improve this procedure would 
consist of trying to find an iterative formula, by which, for eachi,(i=s+1,s+2,...,c), 


Wi would be determined as a function of Wie y> the number of operations required to compute 
{W,} with such a formula would be proportional to c. 
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Unfortunately this scheme is wrong: Elementary examples show that Wi depends 
essentially upon variables other than W;_,. In order to design a more efficient method, it is 
necessary to analyze carefully the set of entities which Wi depends upon. 


A NEW ITERATIVE METHOD 
Let 


REMARK: If the sequence of nonnegative real numbers P.. is nonincreasing and h is 
a real number, then the sequence, 


R, = max (0; P,. - h), 


is also nonincreasing. In fact, either all (P,, - h) are positive (in which case Ry, = P. -h for 
all k), or there exists an integer k' such that P. - h<= 0 if and only if k = k'; then Ry, =0 
for k =k’ and R, = P.. -h for k< k). In either case, the sequence Ri. is nonincreasing. 


Let now {Hi}; (m = 1, 2,..., 8) denote a nonincreasing s-tuple of real numbers Ha 
and {™,,}, (m= 1, 2,..., 8 +1) a nonincreasing (s + 1)-tuple of real numbers Mas Then, 


if R is a real number, the insertion of R into the sequence Le which is represented by the 


mapping of ( 5 H U R onto st M__ and defined by the equations, 
m=1 ™ m=1 ™ 
Mn = Hy whenever m= s and Ti > 8, 
(3.1) M., = Hm-1> Whenever m > 1 and HL, = R, 
Mn = R, in any other case, 


requires an average number of operations proportional to s. 
Let us now go back to Pollaczek's procedure. Let {B AN (m= 1, 2,..., 8) denote the 


sequence consisting of the first s elements of the sequence {Bi}. Hence, 


. gi 
(3.2) W; = B s° 
The analysis of Steps (a), (b) and (c) of the second section herein (Pollaczek's Proce- 
dure and Its Handicaps) shows that {Bit} may be computed from the sequence {B a by 
performing the following two steps: 
‘ , = ba , =e i — wo 
I. Calculating the sequence, {#,, } = max (0; Bn - Aint + A;) = max (0; in 144) 


(m = 1, 2, ..., 8). 
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Il. Inserting the number, 


G;,, = max (0; W, + T; - 144) ; 
into {#i,,} in a way such that the nonincreasing sequencing is preserved. Thusa 
sequence {™,, }(m = 1, 2,..., +1) is obtained. 


It readily follows from the above discussion that the equations, 


(3.3) W.., = M, (i= 8s, s+1,...,c-1), 


hold. This yields the following 


THEOREM: Ww, of is fully determined when the following items are known: 

(1) The s numbers =. (m = 1, 2,..., 8), the last of which is W,. 

(2) The time T; the i-th customer is held for service. 

(3) The time I 41 lapsing between the arrivals of the i-th and of the (i +1)-st 


customer, 


THE SEQUENCE OF OPERATIONS 
Figure 1 shows the flow chart of a subroutine able to perform the operations which are 
necessary for finding the sequence {w,} using the recurrence relations (3.3). 


In Figure 1, Dm,(m = 1, 2, ... , s) denotes the location where Bi is stored, (D,) 
denotes the contents of the location Dw L denotes a working storage location, and (L) 
denotes its contents. 


The procedure, which consists of iterating Steps I and II of the preceding section (A 
New Iterative Method), enables us to compute the whole sequence {W;}, (i= 1, 2,...,¢) with 
a number Po ( c) of operations, such that 


og(c) = K,-c°s, 


where Ky is a positive constant. 
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START ° 
Dm = 0, (m=1 «+, s) 


Wm = 0,(m=1---, $) Dm = max (0; (Dm)-1i ) 
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Figure 1 - Sequence of operations 
































ANNOUNCEMENT 

A Fortran program to compute Negative Binomial Probability Distribution Tables is 
now available. 

The program operates in two modes: (a) input of a single mean and variance to output 
a single table and (b) input of a single mean and initial, incremental, and final variances to 
output a series of tables. Unlimited input in a chosen mode is restricted to one card (record) 
per input. 

Output tables include a heading, page number, associated mean and variance, the 


numerical value of the variables Q, k, 1- = Pi P_, and four successive probabilities across 
j 
each line. Each supplementary page of a table has its own page number and a modified head- 


ing. Up to 1000 probabilities can be computed and printed in floating point mode, from P A to 
Pogg> with a numerical threshold chosen by the user. 

Copies of the program, with operating instructions, will be provided upon request. 

All inquiries should be addressed to: 


oO’ 


Logistics Research Project 

The George Washington University 
707 22nd Street, N. W. 
Washington 7, D. C. 


Attn: M. Hershkowitz 








RECENT PUBLICATIONS 


TESTING STATISTICAL HYPOTHESES, by E. L. Lehmann. John Wiley and Sons, Inc., 
New York, 1959. 


The publication of this book makes generally available Professor Lehmann's excellent 
set of lecture notes on Testing Statistical Hypotheses. The mimeographed notes have proved 
exceedingly useful to those of us fortunate enough to have had them for the past decade. The 
pook is, however, more than just a printed version of the mimeographed notes. There is much 
material and many problems that do not appear in the notes. The rewriting and adding of mate- 
rial has in no way lessened the clarity of Professor Lehmann's exposition, if anything it has 
enhanced it. 

The main purpose of the book is a systematic treatment of the Neyman-Pearson for- 
mulation of the theory of hypothesis testing and confidence sets. This the author does clearly 
and rigorously in the mathematical framework of measure theory in abstract spaces. The 
clarity is indeed a major contribution. It is possible, though difficult, to read much of the text 
with a background that only includes the calculus. The author's chapter 2, 'The Probability 
Background," is brief (29 pages) but reasonably complete as a background for the text. The 
reviewer has considerable doubt that any of the short background chapter in Mathematical 
Statistics Texts do more than set the tone of the text and refresh the trained reader's memory 
slightly. 

There are clearly worked out theoretical illustrations of the subject matter in each 
chapter. There are numerous, well selected, theoretical problems which further illustrate 
and extend the results of each chapter. The author apologizes that the bibliography at the end 
of each chapter is not aimed at completeness. The reviewer only wishes that other authors 
were half as complete. The sections on nonparametric tests are, in the reviewer's opinion, 
some of the clearest treatments ever given to these topics. 

The only complaints the reviewer has in connection with the book are (a) the total lack 
of numerical examples, and (b) the furtherance of the misleading nomenclature that has grown 
up in statistics, e.g., admissibility, most stringent test, etc. Since the author explained his 
view on the first, . . . readers will usually have had previous experience with statistical 
methods, ... ,'' this implies the reviewer has some doubts about the general reader's ability 
in recognizing the applicability of the techniques. On the second point, the author is not, of 
course, responsible for the connotations of some statistical nomenclature or for this reviewer's 
great aversion to the nomenclature. 

The reviewer is, as one can readily gather, enthusiastic about Professor Lehmann's 
book and feels that anyone interested in the theoretical aspects of Statistical Inference will 
benefit from reading it. 


(Reviewed by Dr. William R. Allen) 





RECENT PUBLICATIONS 


MATHEMATICAL METHODS FOR DIGITAL COMPUTERS, by Anthony Ralston and 
Herbert S. Wilf. John Wiley and Sons, New York, 1960, xi + 293 pp. 


This book deals with those topics in numerical analysis which have been influenced 
most by the development of the modern electronic computer. Consequently, the emphasis is 
on new methods which look promising and on newly developed techniques. The book is divided 
into six sections, and it contains 26 chapters in all, each written by an expert in the field. 

The first chapter contains a good description of the various aspects one has to keep in 
mind in searching for the most economical mathematical procedure. There one finds also the 
observation that more recently rational approximations are preferred to polynomial approxi- 
mations. Then the best methods for generating elementary functions are discussed. 

The second section of the book contains several good chapters on computing the inverse 
of a matrix: the Gauss-Jordan method (using the product form of the inverse), the Gauss- 
Seidel method, the conjugate gradient method (diagonalizing a matrix by orthogonal transfor- 
mation), and the method of rank annihilation. There is also a good chapter on the determina- 
tion of the eigenvalues and eigenvectors of a matrix by Jacobi's diagonalization method. 

The third section presents a good description of numerical integration methods for the 
solution of ordinary differential equations; in particular, Milne's iterative method is compared 
with Hamming's noniterative method. Also there is a chapter on the Runge-Kutta-Gil method 
of integration in which differential operators are extensively used but which is, on the whole, 
rather pragmatic. This is followed by a chapter on the numerical solution of boundary value 
problems and one on differential equations with large time constants. 

The fourth section contains an excellent chapter on the solution of parabolic partial 
differential equations by the method of characteristics and two good chapters on the solution of 
hyperbolic partial differential equations. There are also two chapters on the solution of elliptic 
partial differential equations: one by iterative methods, and one by Monte Carlo methods. 

The fifth section deals with statistical topics. It contains a good description of factor 
analysis and good descriptions of autocorrelation and spectral analysis. There are also 
chapters on multiple regression analysis and the analysis of variance. 

The sixth section deals with miscellaneous subjects. It contains a good chapter on the 
numerical solution of polynomial equations by an improved Bernoulli method and a good 
chapter on an improved Gaussian method for numerical quadrature. There are also chapters 
on quadratures by Monte Carlo methods, Fourier analysis, linear programming, and network 
analysis. 

In the beginning of each chapter, the purpose of the computation is indicated, and this 
is followed by a mathematical discussion. In some of the discussions, the theorems are neatly 
stated and elegantly proved; most of them contain also a detailed erroranalysis, as it should be. 
Next the calculation procedure is outlined, and the flow chart is shown and described. This is 
followed by indicating the subroutines that might be required, working out a sample problem, 
stating the memory requirements, and estimating the running time. The list of references at 
the end of each chapter is often extensive, and most of the referenced books, articles, and the 
like, have been published in the last 10 years. In summary, the book probably will be indis- 
pensable to everyone working in the numerical analysis field. Therefore, the book is highly 


recommended to all those interested in digital computers as powerful tools in the numerical 
solution of problems. 


(Reviewed by Herman F. Karreman) 
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PROBABILITY, AN INTRODUCTION, by Samuel Goldberg, Prentice-Hall, Inc., 1960, 
ix + 322 pp., $7.95. 


"This book is intended for all who require a mathematically sound but elementary 
introduction to the theory of probability." 

With this statement of intent, the author begins a very carefully written book. In the 
preface, it is further explained that the author's purpose is to provide a text for a one-quarter 
course in probability for those who do not have a working knowledge of calculus. To do this, 
the choice of topics has been limited to those involving only finite sample spaces. A brief 
outline of the book's five chapters follows. 

In Chapter 1, the author carefully introduces the basic concepts of elementary set 
theory. The algebra of sets is discussed quite formally, with emphasis on membership (truth) 
tables and Venn diagrams. The last section studies ordered pairs and Cartesian products. 

Finite probability spaces are studied in Chapter 2. After defining and illustrating the 
notions of sample space and event, a probability distribution is defined formally as an assign- 
ment of nonnegative numbers to the simple events such that the sum of these numbers is one. 
The principle properties of probability distributions are then proven, followed by a study of 
conditional probabilities, and a detailed discussion of independence. The chapter concludes 
with a probability model in genetics. 

A short Chapter 3, entitled Sophisticated Counting, follows in which are derived the 
number of r-tuples from n objects, the number of subsets of size r from n objects and the 
number of distributions of n similar objects among k boxes. Logarithms are used as an aid 
in working examples. The Binomial and Multinomial expansion theorems are proven and 
applied. 

Random variables are introduced, in Chapter 4, as functions having a sample space for 
domain and a set of real numbers for range. Distribution function, mean and variance are 
defined and the linear properties of expectations derived. Joint distributions, covariance and 
correlation are discussed. Chebyshev's inequality and the law of large numbers are proven. 
The fact that a correlation coefficient is + 1 if and only if the two random variables concerned 
are linearly related is also proven. 

The last chapter is about Bernoulli trials and the Binomial distribution. The mean and 
variance of the Binomial distribution are derived. A two-page table of cumulative Binomial 
probabilities is provided. Hypothesis testing is discussed in detail in this Binomial context. 

In the preface the author suggests three possible audiences which might profitably use 
his book. They are (a) high school teachers without a calculus prerequisite but who may be 
called upon to teach a twelfth grade course in probability and statistical inference such as that 
recommended by the Commission on Mathematics of the College Entrance Examination Board; 
(b) college students who are required to take a one-year course in mathematical thinking, the 
author believing that ''a strong case can be made for a course that concentrates on sets and 
probability in the finite case at first, proceeds to an introduction to the calculus, and then 
applies this calculus to the elements of probability in the infinite case"'; (c) college students 
with calculus who might use this text in an introductory course in probability, designed as a 
motivational preliminary to the usual mathematical statistics course. Concerning the sentence 
quoted at the outset of this review, it may well be asked whether the set of people who "require 
both a "mathematically sound" and an "elementary" course in probability is non-empty. That 
it is non-empty is evidenced best by the existence of audience (b) just mentioned, since many 
universities do require of a large proportion of the freshmen to take a course in modern 
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fundamental mathematics, and since probability theory does present an easily motivated frame- 
work within which to present rigorous deductive reasoning. There is some question in the 
reviewer's mind about the sizes and importance of audiences (a) and (c). Concerning (a), it 
would be desirable and should be encouraged that teachers of a precalculus course in proba- 
bility have a calculus prerequisite. There is also the question as to whether or not a twelfth 
grade course which would almost always be a terminal course should emphasize the inductive 
reasoning of statistical inference rather than the careful deductive reasoning of Goldberg's 
book. Nevertheless, subject to the intent of the author and the suggested audiences, it must be 
said that the book has been very carefully and clearly written. Errors and misprints are 
exceptionally rare. This is undoubtedly due in great part to the fact that preliminary forms of 
this manuscript were used in several courses at Oberlin College for 2 or 3 years prior to 
publication which is a very desirable but often neglected procedure in preparing a book. 

The book contains 360 problems together with solutions to half of them. Many of the 
problems, particularly in Chapters 1 and 3, are very interesting. This reviewer regrets, 
however, the author's policy of introducing important concepts in the exercises. Approxi- 
mately 22 concepts of different degrees of importance are so introduced, including maximum 
likelihood estimation, operating characteristic curve (which is not linked up with the power 
function introduced elsewhere), and the hypergeometric distribution. 

In the reviewer's opinion, there are very few criticisms that can be made of the author's 
presentation. This reviewer would have preferred having probabilities for equally likely events 
emphasized and introduced prior to the general definition on page 55. The author considers the 
equally likely case much later on page 69, although 7 of the 9 examples occurring between pages 
55 and 69 concern equally likely probability models, as do 15 out of the 17 numerical exercises 
related to this portion. On page 70, the author emphasizes that "'this rule for computing prob- 
abilities (namely, number of outcomes in an event over the total number of outcomes) is appli- 
cable only when all simple events have been assigned the same probability . . . and does not 
apply to the wide variety of important problems where it is not reasonable to make this special 
assignment."" However, it seems impossible for the student reader to ascertain and acknowl- 
edge this ''wide variety of important problems," It is also felt that the main properties of 
probabilities should be motivated by the analogous properties of frequencies, the latter not 
being mentioned at all in this text. The author uses only the undefined term "natural" to 
describe many equally likely probabilities (e.g., those for fair coins, fair dice, random card. 
selecting and random number drawing). In an occasional place, concepts are used before they 
are defined. For example, "equally likely" is defined on page 69 but used as early as pages 
56 and 59, a "fair" coin is defined on page 61 but used on pages 56 and 60, independent random 
variables are defined on page 204 but referred to on page 202. Chebyshev's inequality is more 
instructively and simply proved as a corollary of the result given as Exercise 3.20 on page 197. 
This reviewer also prefers to have random variables studied simultaneously with probabilities, 
rather than confined to a later chapter as in this text. 

In summary, this book should be given very serious consideration by any instructor 
giving a mathematically sound, elementary course in finite-sample-space probability theory. 


(Reviewed by Ronald Pyke) 





frame- 
he 

, it 
Oba- 
elfth 
ictive 
z's 
ust be 


> 


‘ms of 


juthor's 
events 
ors the 
Nn pages 
rcises 
prob- 
appli- 
not 
pecial 
owl- 
f 
ot 


ard. 
they 
eS 
ndom 
more 
ze 197, 
lities, 


RECENT PUBLICATIONS 


HANDBOOK OF AUTOMATION, COMPUTATION, AND CONTROL, by Eugene M. 
Grabbe, Simon Ramo, and Dean E. Woolridge. John Wiley and Sons, Vol. I, 1958, xx + 1020pp., 
Vol. II, 1959, xxii + 1093 pp. 


The rapid development of automatic control systems during and after World War II 
created a vast body of technical knowledge in this area. Because of its complexity, there 
arose the need for codification and unification of this knowledge. The two-volume handbook 
fills this need and provides an excellent guide to work that has been done; it gives the systems 
engineer the basic information on new techniques and components that he needs for designing 
and developing control systems in general. 

Volume I, consisting of five parts, is devoted to Control Fundamentals. Part A contains 
a selection of mathematical subjects which are essential for a good understanding of what fol- 
lows. Some of these subjects—e.g., sets and relations, Boolean algebra, probability, and 
statistics—are not ordinarily found in engineering handbooks. Part B covers, to some extent, 
the numerical analysis with which all those interested in the efficiency and accuracy of com- 
puting performances should become familiar. Part C deals with various aspects of operations 
research, a subject that lies, in this reviewer's opinion, somewhat outside the scope of this 
Handbook. Part D gives an up-to-date account of information theory, smoothing and filtering 
of data, and data transmission. Part E, occupying the second half of Volume I, covers rather 
thoroughly the present stage of feedback control. It consists of eight chapters. The first two 
chapters deal with the methodology of feedback control and the fundamentals of system analysis. 
The next several chapters deal with stability, relationships between transient and frequency 
response, noise and random inputs, etc. This part (Part E) and Part D are undoubtedly the 
most interesting ones, viewed from the standpoint of the systems engineer. 

Volume II, also consisting of five parts, is entirely devoted to the electronic computer 
as a data processing device, and it occupies a central position in this Handbook. An intro- 
ductory part on computer terminology is followed by three parts on the digital computer. In 
Part B, which deals with programming, the subject of automatic programming is well covered 
in Sections eight to 15 of Chapter 2. In these sections, one finds a detailed description of 
assembly programs, sunroutine and utility programs, and compilers and translators. A Soviet 
algebraic language compiler is discussed in Section 13; this is followed by a section on inter- 
preters and one on recursive languages. Part C is concerned with the use of the digital com- 
puter for purposes of data processing. After a discussion of the requirements that have to be 
met, examples of applications in business (Chapters 8 and 9), as well as science and engineer- 
ing (Chapter 10), are presented. Next there is a chapter on the use of the digital computer to 
perform logical operations and to translate languages. Part D discusses in great detail the 
design of the digital computer and is, of course, again of particular interest to the systems 
engineer, After an introductory chapter on design fundamentals, various techniques for 
increasing the reliability of the computer are discussed in Chapters 13 and 14. This material 
is followed by a chapter on magnetic cores and one on translators—the latest developments in 
this area. Chapter 17 deals with the logical design of the digital computer and contains a brief 
discussion of the Boolean algebra on which it is founded. Section 5 of this chapter shows how a 
computer can be used for the detection of mistakes in the logical design of a new computer by 
way of simulation. Part E is devoted to the other main type of electronic computer—the analog 
computer; it contains a chapter (Chapter 24) on the use of this computer in the study of physical 
Systems and one (Chapter 25) on the solution of engineering problems. Chapter 26 discusses 
Statistical applications and random-noise theory. Finally, Part F deals with some unusual 
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computer systems, such as combined analog-digital computers and simple turing-type com- 


puters. 

There is an extensive list of references at the end of each chapter for the reader who 
wishes to delve deeper and learn more on the subject of his particular interest. This makes 
owning this Handbook all the more worthwhile for all those genuinely interested in the related 
subjects of automation, computation, and control. 


(Reviewed by Herman F. Karreman) 


* U.S. GOVERNMENT PRINTING OFFICE : 1961 O—609176 
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